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1.  EXECUTIVE  SUMMARY 


This  report  describes  work  accomplished  under  contract  option  1  of  the  EW  Receiver  and 
Processing  Concepts  Evaluation  Program  (RAPCEval  2).  The  Air  Force  Research  Laboratory 
(AFRL)  awarded  this  contract.  This  branch  of  AFRL  is  located  at  Wright  Patterson  Air  Force 
Base,  in  Dayton,  Ohio.  The  Option  1  task  of  this  program  was  awarded  to  Mercer  Engineering 
Research  Center  (MERC)  on  October  15, 1997.  Work  was  completed  on  the  RAPCEval  1  contract 
in  October,  1996,  and  this  contract  was  initiated  in  November  1996.  This  report  describes  work 
accomplished  on  the  contract  from  December  2,  1997  through  March  31,  1999.  The  project 
includes  five  options  and  has  maximum  duration  of  60  months  from  the  date  of  award.  In  support 
of  the  contract,  a  number  of  activities  and  projects  have  been  initiated.  Some  of  these  are  complete, 
and  many  more  are  continuing.  A  number  of  activities  have  spawned  new  directions  for 
investigative  effort. 

The  Program  Research  Standards  Committee  (PRSComm),  established  at  the  outset  of  the  Basic 
task,  has  continued  to  meet  regularly  and  has  provided  valuable  guidance  and  suggestions  both  as 
to  the  direction  of  students  involved  in  the  effort,  and  the  scope  and  emphasis  of  the  research  efforts 
in  general.  Current  membership  of  this  committee  represents  the  Air  Force  Research  Laboratory, 
Robins  Air  Force  Base,  Mercer  University,  and  Mercer  Engineering  Research  Center.  The  specific 
members  are  listed  in  another  section  of  this  report  (section  4,  Research  Support). 

Three  meetings  of  the  PRSComm  were  held  during  contract  Option  1.  Eight  student  research 
reviews  were  presented  and  were  approved  by  the  PRSComm  in  die  course  of  these  meetings  (see 
section  5,  Project  Activity).  Valuable  discussion  and  suggestions  were  provided  for  directing  and 
focusing  the  students'  work.  In  addition,  suggestions  were  received  for  topics’  reduction  in  scope 
or  better  focus.  Frequently,  unanticipated  or  unknown  resources  and  techniques  were  pointed  out 
for  the  benefit  of  the  students.  One  meeting  involved  a  presentation  by  AFRL  researcher  James 
Stephens,  who  delivered  a  visiting  lecture  with  an  emphasis  on  communication  countermeasures, 
with  focus  on  advanced  digital  signal  processing  for  electronic  warfare. 

Students  engaged  in  the  RAPCEval  program  have  continually  enjoyed  fruitful  contact  with 
knowledgeable  personnel  at  the  Air  Force  Research  Laboratory  in  their  respective  areas  of  interest. 
They  have  also  interacted  with  experienced  colleagues  at  Robins  Air  Force  Base,  employees  of 
MERC,  Mercer  University  School  of  Engineering  faculty,  and  various  representatives  of  industry. 

It  should  be  noted  that  software  generated  under  this  contract  is  not  government-owned. 

The  RAPCEval  contract  has  stimulated  gratifying  communication  and  collaborative  research 
effort  among  students,  university  faculty,  MERC  personnel,  personnel  at  the  Air  Force  Research 
Laboratory  and  Warner  Robins  AFB  and  industry.  All  parties  have  expressed  satisfaction  with 
the  contract  results. 
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2.  INTRODUCTION 

These  tasks  specify  requirements  for  analytical  and  research  support  of  in-house  research  at  the 
Air  Force  Research  Laboratory,  AFRL/SNR.  There  is  increasing  sophistication,  quantity,  and 
mobility  of  hostile  radars,  AAM,  SAM,  and  AAA  fire  control  systems.  Electronic  warfare 
receivers  for  radar,  electro-optic,  infrared,  ultraviolet  missile  warning  and  electronic 
countermeasures  need  operational  upgrades  to  allow  penetrating  aircraft  acceptable  survivability. 
This  encourages  maintenance  of  in-house  labs  to  support  development,  to  evaluate  concepts,  and 
to  test  new  receivers,  processors,  and  software. 

2.1  EW  Receiver  Effort 

Complex  EW  environments  have  caused  employment  of  numerous  receiving  systems. 
Augmentation  of  in-house  capability  for  evaluation,  novel  concept  development,  and  exploitation 
of  new  technology  is  needed.  Computer-aided  simulation  of  new  systems  and  concepts  can  save 
resources.  New  high  speed  analog-to-digital  (A/D)  converter  technology  may  allow  input 
frequencies  to  be  digitized  in  baseband  before  the  crystal  video  detector,  possibly  allowing  real 
time  digitized  frequency,  pulse  width,  and  pulse  amplitude.  Advancing  materials  technology  for 
infrared  (IR)  /  ultraviolet  (UV)  /  radio  frequency  (RF)  energy  offers  the  possibility  of  augmented 
and  combined  sensors.  Investigation  of  these  materials  is  needed  to  reduce  the  kind  and  number 
of  avionics  needed  in  combat.  Filtering  and  discrimination  advances,  both  in  hardware  and 
software,  may  allow  enhancement  of  fielded  EW  systems. 

2.2  EW  Processing  Effort 

A  modem  EW  system  must  face  an  increasing  number  of  hostile  threats  that  are  multimode. 
Sensors  to  intercept  such  threats  now  include  radar  warning  and  electronic  intelligence  systems. 
Information  from  these  sensors  must  be  processed,  the  threat  identified,  and  appropriate 
countermeasures  initiated  to  counter  these  threats.  An  augmentation  of  in-house  capability  is 
required  to  evaluate  processor  hardware  and  software,  to  exploit  novel  ideas,  and  to  investigate 
advanced  concepts  such  as  artificial  intelligence  to  determine  the  nature  of  the  threat,  and  what 
countermeasures,  if  any,  to  employ. 

2.3  EW  Exciter  Effort 

Digital  exciters  are  being  developed  to  provide  a  flexible  active  ECM  asset  against  a  wide  variety 
of  modem  threats.  The  need  exists  to  evaluate  the  various  exciter  architectures,  advance  and 
develop  unique  concepts,  and  advance  the  digital  exciter  technology  base.  The  novel  concepts 
and  technologies  must  be  evaluated  for  effectiveness  against  the  proposed  application. 

2.4  EW  Antenna  Effort 

The  role  of  antennas  as  the  "eyes  and  ears"  of  the  sensor  suites  continues  to  make  RF  antenna 
technology  development  vital  to  the  Air  Force  mission.  Airborne  antenna  apertures  of  the  future 
will  be  low  cost,  broadband,  low  radar  cross  section  (RCS),  and  multifunction  in  nature  to  earn 
their  way  onto  platforms  where  space  is  at  a  premium. 
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3.  SCOPE 

The  overall  program  consists  of  a  basic  task  and  four  options  that  are  the  conglomerate  of 

different  work  efforts  and  technologies  within  the  Electronic  Warfare  arena.  Detailed 

descriptions  are  given  as  follows: 

•  Basic  Task  -  The  basic  task  will  provide  the  tasks  necessary  to  analyze  software  and 
hardware  approaches  to  perform  the  exploratory  development  of  EW  technology  in  these 
technology  areas:  radar  hardware,  laser  hardware,  infrared  hardware,  and  ultraviolet 
hardware.  The  task  will  analyze  receiver  and  exciter  technology  to  generate  ECM  signals  to 
improve  ECM  system  performance.  In  addition,  the  scope  of  the  basic  task  will  include 
signal  processing  technology  related  to  the  hardware. 

•  Option  1  -  These  tasks  will  be  those  necessary  to  analyze  receiver  technology  for  application 
to  modem  digital  spectrum  estimation  techniques  in  order  to  improve  EW  /  SIGINT  /  ELINT 
/  IR  /  EO  receiver  performance. 

•  Option  2  -  This  option  consists  of  those  tasks  necessary  to  identify  high  risk  design  areas  for 
an  EW  /  SIGINT  /  ELINT  /  IR  /  EO  hardware  approach,  to  perform  exploratory  design 
assessments  for  selected  functions,  and  to  determine  the  degree  of  parallel  processing 
achievable. 

•  Option  3  —  This  option  is  “reserved.” 

•  Option  4  -  These  tasks  are  those  essential  to  EW  /  SIGINT  /  ELINT  /  IR  /  EO  hardware  and 
signal  processing  including  but  not  limited  to  pulse-deinterleaving,  parametric  extraction,  and 
threat  identification. 
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4.  RESEARCH  SUPPORT 


For  support  of  the  overall  contract,  a  "Program  Research  Standards  Committee"  has  been 
established.  Membership  for  this  committee  was  most  recently  updated  March,  1 997.  Current 
members  are: 

•  from  the  Air  Force  Research  Laboratory  in  Dayton, 

•  Mr.  Nicholas  Pequignot  (the  program  manager  for  AFRL) 

•  Mr.  Emil  R.  Martinsek 

•  Mr.  Norman  A.  Toto 

•  Dr.  Duane  A.  Warner 

•  Mr.  Paul  J.  Westcott 

•  from  Warner  Robins  Air  Logistics  Center, 

•  Mr.  Steve  Strawn  (the  program  manager  for  WR-ALC) 

•  Mr.  John  LaVecchia 

•  Mr.  Phil  Oliver 

•  Mr.  Ches  Rehburg 

•  Mr.  Larry  Sheets 

•  from  Mercer  University  and  MERC, 

•  Dr.  Tom  Bass  (the  program  manager  for  MERC) 

•  Dr.  David  Barwick,  (chairman  of  the  standards  committee) 

•  Dr.  Aaron  Collins  (Mercer  U) 

•  Dr.  Behnam  Kamali  (Mercer  U) 

•  Dr.  Paul  MacNeil  (Mercer  U) 

The  EW  Receiver  and  Processing  Concepts  Evaluation  Program  was  awarded  to  the  Mercer 
Engineering  Research  Center  (MERC)  by  WPAFB/WL  under  contract  F09603-93-G-0012-0017. 
This  contract  is  administered  through  WR-ALC.  The  overall  program  has  a  funding  ceiling  of 
$499,940.  Incremental  funding  will  be  accomplished  via  a  series  of  contract  options.  The  basic 
contract  is  $99,998.00,  option  (1)  is  $99,998.00,  option  (2)  is  $99,998.00,  option  (3)  is 
“reserved”,  and  option  (4)  is  $50,000.00. 

Contracts  have  been  awarded  for  the  Basic  program,  Option  1,  Option  2,  and  Option  4. 
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5.  PROJECT  ACTIVITY 


5.1  Steering  Committee,  April  1998 

5.1.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 
_ Minutes  2  April  1998 _ 


A  meeting  of  the  Program  Research  Standards  Committee  (PRSC)  for  the  EW  Receiver 
and  Processing  Concepts  Evaluation  Program  (RAPCEval2)  was  hosted  by  Mercer  Engineering 
Research  Center  at  8:30  a.m.  in  the  new  MERC  auditorium.  Committee  members  present  were 
Dave  Barwick,  Tom  Bass,  Aaron  Collins,  Benham  Kamali,  Nick  Pequignot,  Norman  Toto, 
Duane  Warner,  Paul  MacNeil,  Phil  Oliver,  and  Ches  Rehberg.  Also  present  were  several 
representatives  from  Mercer  University,  students  scheduled  to  speak,  and  other  students 
considering  participation  in  the  RAPCEval  program.  Five  student  presentations  were  given. 

Four  were  reports  on  research  in  progress,  while  one  was  a  research  proposal. 

After  welcoming  remarks  by  Dr.  Barwick,  Dr.  Tom  Bass  gave  a  brief  description  of  the 
RAPCEval  program  and  introduced  the  students  who  were  scheduled  to  speak.  The  students  and 
their  topics  are  listed  in  the  table. 


Ron  Brinkley 

“Burst  Error  Correction  with  Reed- 
Solomon  Codes” 

Dennis  Ludwig 

“Identification  of  IR  Images  Using  Neural 
Networks” 

Henderson  Benjamin 

“Selection  of  Reed  Solomon  Codes  Using 
Neural  Networks” 

Randy  Ford 

“Passive  Location  via  Evolutionary 

Genetic  Algorithms” 

Dennis  Ludwig 

“Identification  of  IR  Images  Using  Neural 
Networks” 

“Hardware  Exploitation  of  RAD  Filtration” 

Ron  Brinkley  showed  overheads  introducing  the  concept  of  Reed-Soloman  error 
correction  methods,  pointing  out  that  RS  codes  are  block  codes.  He  also  gave  a  summary  of  the 
types  of  communication  errors,  together  with  interleaving  techniques  designed  to  combat  the 
errors.  Ron  discussed  his  analysis  and  tabulation  of  the  Reed-Soloman  error-correcting  codes  for 
specific  classes  and  burst  error  categories.  Announcement  was  made  of  Ron’s  presentation  on 
this  subject  to  the  IEEE  Pacific  Rim  Conference  (PACRIM  ’97). 
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Henderson  Benjamin  reviewed  his  topic  on  the  application  of  neural  networks  to  the 
selection  of  Reed-Solomon  (RS)  codes.  Henderson  showed  many  of  the  details  of  the  manner  in 
which  RS  codes  are  selected  by  a  neural  net  (NN).  Included  were  training  techniques  for  the  NN 
and  performance  statistics  for  the  trained  network.  Henderson  acknowledged  dependence  on  the 
work  of  Ron  Brinkley  for  availability  of  code  selection  tables  that  were  required  for  the  NN 
application.  This  work  is  nearly  complete  and  the  thesis  is  under  final  review.  An  abbreviated 
version  of  the  thesis  is  available  in  the  RAPCEval  annual  report,  dated  3  February  1998. 

Dennis  Ludwig  reviewed  the  pose  estimation  problem  for  an  object  identified  by  an 
infrared  (or  other)  image.  He  had  planned  to  extend  a  two-dimensional  solution  developed  at  the 
Wright  Laboratory  to  a  general  solution  in  three  dimensions.  The  chief  difficulty  encountered  is 
the  necessity  of  taking  available  two-dimensional  images,  and  somehow  recreating  the  third 
dimension.  A  number  of  questions  were  evoked  by  this  presentation,  ranging  from  whether 
three-dimensional  data  is  available,  to  the  specifics  of  how  the  images  are  brought  into 
alignment. 

Tracy  Tillman  presented  additional  details  on  the  development  of  a  unique  method  of 
table  generation  and  table  lookup  for  use  in  RAD  parameter  production.  He  has  discovered  some 
limitations  for  the  technique,  and  plans  to  finish  a  report  on  his  discoveries,  calculations,  and 
algorithmic  development  in  the  next  quarter  or  two. 

Dr.  Barwick  chaired  the  committee  discussions  that  followed  the  last  presentation. 
Suggestions  were  gathered  regarding  the  student  presentations.  The  group  agreed  that  the 
communication  research  projects  directed  by  Dr.  Kamali  were  very  pertinent  to  current  AF  needs 
and  should  be  continued.  Several  sensed  that  Dennis  Ludwig  seemed  to  have  lost  focus  of  the 
original  goals,  probably  due  to  the  departure  of  his  thesis  advisor  recently.  Dr.  Bass,  in 
conjunction  with  Dr.  Collins,  will  be  giving  increased  direction  and  steer  him  back  on  track.  A 
question  arose  as  to  whether  the  RAD  lookup  tables  were  in  use  in  the  PLAID  program. 
Currently,  they  are  not,  but  further  evaluation  is  planned  during  the  next  phase  of  R&D.  The 
comment  was  made  that  Randy  Ford’s  schedule  seemed  optimistic.  Another  question  concerned 
the  potential  to  test  the  differential  evolution  solution  technique  on  actual  PLAID  data.  This  will 
be  investigated. 

During  final  discussions,  Norm  Toto  said  that  topics  formerly  being  investigated  at  AFIT 
were  “going  away.”  He  thought  MERC  might  be  able  to  continue  some  of  these  via  RAPCEval. 

Norm  also  suggested  potential  RAPCEval  research  investigations  into  GPS  anti-jamming 
techniques  and  related  technologies. 

At  this  point,  the  meeting  was  adjourned. 
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5.1.2  Meeting  Agenda 


Agenda  steerin 

m 

RAPCEval 

nrnnr^m  Main  Auditoriur 

program  135( 

G  COMMITTEE 
EETING 

98  April  2  —  8:30  AM  to  2:30  PM 
n,  Mercer  Engineering  Research  Center 
Osigian  Boulevard,  Warner  Robins,  GA 

Meeting  till;  Nicholas  Pequignot,  AFRL/SNRP,  Facilitator:  Dr.  Tom  Bass 

called  by:  ;  Air  Force  Research  Laboratory  Program 
y  Manager 

Committee  AF  Rsch  Laboratory 

Members:  Mr  Nicholas  Pequignot 

Mr  Emil  R  Martinsek 

Mr  Norman  A  Toto 

Dr  Duane  A  Warner 

Mr  Paul  J  Westcott 

_  .j  1  pequigna@sensors.wpafb.af.  mil 

t  m  3*1  '•  martinser@sensors.wpafb.af.  mil 

addresses:  totona@sensors.wpafb.af.mil 

wamerda@sensors.wpafb.af.mil 
::•••  westcopj@sensors.wpafb.af.mil 

Robins  AFB 

Mr  Steve  Strawn 

Mr  John  LaVecchia 

Mr  Phil  Oliver 

Mr  Ches  Rehberg 

Mr.  Larry  Sheets 

lnerr@ewir-wr.robins.af.mil 

lavecchia@ewir-wr.robins.af.mil 

oliver@ewir-wr.robins.af.mil 

rehberg@ewir-wr.robins.af.mil 

sheets@ewir-wr.robins.af.mil 

Mercer  University 

Dr  David  Barwick 

Dr  Tom  Bass 

Dr  Aaron  Collins 

Dr  Benham  Kamali 

Dr  Paul  MacNeil 

dbarwick@merc.mercer.edu 
tbass@merc.mercer.edu 
collins_as@mercer.edu 
macneiljpe@mercer.peachnet.edu 
kamal  i_b@mercer.peachnet.edu 
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Meeting  Agenda,  continued. . . 


Schedule 

Continental  Breakfast  /  Greetings 

Dr  David  Barwick 

8:30  AM -9:00  AM 

Meeting  Overview 

Dr  Tom  Bass 

9:00  AM -9:20  AM 

Student  Proposals 

Dr  Tom  Bass  /  Students 

Burst  Error  Correction  with 
Reed-Solomon  Codes 

Ron  Brinkley 

9:20  AM  -9:55  AM 

Selection  of  Reed  Solomon 

Codes  Using  Neural  Networks 

Henderson  Benjamin 

9:55  AM  -  10:30  AM 

Break 

10:30  AM -10:50  AM 

Student  Proposals  (continued) 

Dr  Tom  Bass  /  Students 

Passive  Location  via  Evolution¬ 
ary  Genetic  Algorithms 

Randy  Ford 

10:50  AM- 11:25  AM 

Identification  of  IR  Images  Us¬ 
ing  Neural  Networks 

Mr  Dennis  Ludwig 

11:25  AM -NOON 

Lunch 

NOON  -  1:00  PM 

Student  Proposals  (continued) 

Dr  Tom  Bass  /  Students 

Hardware  Exploitation  of  RAD 
Filtration 

Mr.  Tracy  Tillman 

1:00  PM- 1:35  PM 

Discussions  &  New  Business 

Dr.  Tom  Bass 

1:35  PM -2:30  PM 

Adjourn 

2:30  PM 
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Toto,  Norman  A.  AFRL/SNRW  (937)  255-4933  toto@sensors.wpafb.af.mil 

Warner,  Duane  AFRL/SUJO  (937)255-4174  warnerda@aa.wpafb.af.mil 


5.1 .4  Overview  of  the  Program  (Dr.  Bass) 

The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  7  pages. 
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April  2, 1998 


All  research  has  been  determined  by  the  committe 
to  be  useful  to  the  Air  Force,  as  well  as  having  me 
for  graduate  research 
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Kamali,  Dr.  Paul  MacNeil 


00 
fa  ON 
P-i  Os 
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Tracy  Tillman,  “Hardware  Exploitation  of  RAD 
Filtration” 
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Improved  Communication  of  Time  and  Time 
Measurements”,  unclassified  report  available  from  MERC 
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★  Kirk  Wright,  “Object  Oriented  Modeling  of  the 

AN/ALQ-172”,  publication  in  process 
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Networks 

Tracy  Tillman,  “Hardware  Exploitation  of  RAD  Filtration 


5.1.5  Presentation  By  Ron  Brinkley 


The  student  briefing  presented  by  Ron  Brinkley  at  this  meeting  is  reproduced  on  the  next  32 
pages. 
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Research  Topic:  BURST  ERROR  CORRECTION  WITH  REED-SOLOMON  CODES 
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_  Ronald  L.  Brinkley 

RF  RECEIVER  AND  PROCESSING  MSEE  Program 

CONCEPTS  EVALUATION  PROGRAM  Final  Briefing 
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[Ski]  B.  Sklar,  Digital  Communication,  Fundamentals  and  Applications,  Englewood  Cliffs,  NJ:  Prentice  Hall,  1988. 
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Channel  Memory 
Long  Single  Burst 
Multiple  Short  Burst 


Convolutional  Interleaving 

•  More  Efficient  than  Block  Interleaving 

•  CIRC 


BLOCK  INTERLEAVER 


Interleaver 


Input  Stream:  Al,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  A10,  A1 1,  A12,  Bl,  B2,  B3,  B4,  B5,  B6, 


.,F9,  FI0,F1 1,  F12 


4 -  n  symbols  - ► 


Al 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

A9 

A10 

All 

A12 

Bl 

B2 

B3 

B4 

B5 

B6 

B7 

B8 

B9 

BIO 

Bll 

B12 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

CIO 

Cll 

02 

D1 

D2 

D3 

D4 

D5 

D6 

D7 

D8 

D9 

DIO 

Dll 

D12 

El 

E2 

E3 

E4 

E5 

E6 

E7 

E8 

E9 

E10 

Ell 

E12 

FI 

F2 

F3 

F4 

F5 

F6 

F6 

F8 

F9 

F10 

FI  1 

F12 

Output  Stream: 


mm,  ci,  Bl,  El,  F1,A2,  B2,C2,  D2,  E2/F2 


A3,  B3,  C3,  D3,E3,F3,A4, 


C12,  D12,  E12,  F12 


4> 


burst 


-► 


De-interleaver 


Input  Stream  Al&%  Cl,  D1,E1,  FI '  A2,B2,  C2,  D2,  E2,E2,  A3,  B3,  C3,  D3,  E3,  F3,  A4, 


4 -  n  symbols  - ► 


At 

A2 

A3 

A4 

A5 

A6 

A7 

A8 

A9 

A10 

All 

A12 

Bl 

B2 

B3 

B4 

B5 

B6 

B7 

B8 

B9 

BIO 

Bll 

B12 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

CIO 

Cll 

02 

m 

D2 

D3 

D4 

D5 

D6 

D7 

D8 

D9 

DIO 

Dll 

D12 

El 

E2 

E3 

E4 

E5 

E6 

E7 

E8 

E9 

E10 

Ell 

E12 

FI 

'F2’ 

F3 

F4 

F5 

F6 

F6 

F8 

F9 

F10 

Fll 

F12 

02,  D12,  E12,  F12 


Output  Stream: 


A1,A2 


A3,  A4,  A5,  A6,  A7,  A8,  A9,  A10,  Al  1,  A12, 


B1B2 


B3,  B4,  B5,  B6, . ,F9,  F10,  FI  1,  F12 
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[Fori]  G.  D.  Forney,  Jr.,  “Burst  Correcting  Codes  for  the  Classic  Bursty  Channel,” 

IEEE  Transactions  on  Communication,  Volume  COM-19,  pp.  772-781,  October  1971. 


_ _  _  _  _  _  Ronald  L .  Brinkley 

RF  RECEIVER  AND  PROCESSING  MSEE  Program 

CONCEPTS  EVALUATION  PROGRAM  Final  Briefing 
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RF  RECEIVER  AND  PROCESSING 
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_ _  Ronald  L.  Brinkley 
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[Kam]  B.  Kamali,  Lecture  notes  on  Reed-Solomon  Codes,  EEE  691.80  Spring  Qrt.  1995,  Mercer  University. 
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699125*0 
9605 WO 
Z2989S0 
ZS  1.265*0 
989519*0 
912699*0 
SW.299'0 
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K)860Z*0 

111 
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226908*0 
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mm 

ita-Hiaaai 


(255,251) 


ItaHgHlI 


(255,247) 


(255,245) 


(255,243) 


(255,241 


(255,239) 


(255,237) 


(255,235) 


(255,233) 


(255,231) 


(255,229) 


(255,227) 


(255,225) 


(255,223) 


(255,221) 


(255,219) 


(255,217) 


(255,215) 


255,213 


(255,211) 


(255,209) 


(255,207) 


(255,205) 


(255,203) 


(255,201) 


255,199 


(255,197) 


(255,195) 


(255,193) 


(255,191) 


(255,189) 


(255,187) 


255,185 


(255,183) 


(255,181) 


(255,179) 


(255,177) 


(255,175) 


(255,173) 


255,171 


(255,169) 


(255,167) 


(255,165) 


(255,163) 


(255,161) 


(255,159) 


(255,157 


(255,155) 


(255,153) 


(255,151) 


(255,149) 


(255,147) 


(255,145) 


(255,143) 


(255,141) 


(255,139) 


(255,137) 


(255,135) 


(255,133) 


(255,129) 


0.984314 


0.976471 


0.968627 


0.960784 


0.952941 


0.945098 


0.937255 


0.929412 


0.921569 


0.913725 


0.905882 


0.898039 


0.890196 


0.882353 


0.874510 


0.866667 


0.858824 


0.850980 


0.843137 


0.835294 


0.827451 


0.819608 


0.811765 


0.803922 


.796078 


0.788235 


0.780392 


0.772549 


0.764706 


0.756863 


.749020 


0.741176 


0.733333 


0.725490 


0.717647 


0.709804 


0.701961 


0.694118 


0.686275 


0.678431 


0.670588 


0.662745 


0.654902 


0.647059 


0.639216 


0.631373 


0.623529 


0.615686 


0.607843 


0.600000 


0.592157 


0.584314 


0.576471 


0.568627 


0.560784 


0.552941 


0.545098 


0.537255 


0.529412 


0.521569 


0.513725 


0.505882 


0.562500 


0.708333 


0.781250 


.825000 


0.854167 


0.875000 


0.890625 


0.902778 


0.912500 


0.920455 


0.927083 
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5.1.6  Presentation  by  Henderson  Benjamin 

The  student  briefing  presented  by  Henderson  Benjamin  at  this  meeting  is  reproduced  on  the  next 
24  pages. 
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_ _ _  Henderson  C.  Benjamin 
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Collect/ Analyze  data  for  shortened  RS  codes  of  length  255 
Collect/Analyze  additional  data  deemed  useful  to  a  designer 
Develop  Neural  Network 
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ing,  and  recording  systems. 
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in  interleaved  cases 
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The  trade  off  of  shortened  codes  is  .the  efficiency  of  the  coderate  10  and  12% 
respectively  in  the  above  case. 
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The  current  cycle  time  of  a  neuron  is  approximately  a  millisecond  which  is  a 
slower  than  semiconductor  devices. 
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approximation  of  Newton’s  method. 
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&4.40%  Correct  for  the  code  dimension 


5.1.7  Presentation  by  Randy  Ford 

The  student  briefing  presented  by  Randy  Ford  at  this  meeting  is  reproduced  on  the  next  26  pages. 
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Research  Topic: 

COMPARISION  OF  DIFFERENTIAL  EVOLUTION  TO  THE  SIMPLEX  METHOD  IN 
OPTIMIZATION  DURING  PASSIVE  EMITTER  LOCATION. 
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Differential  Evolution,  developed  by  Kenneth  Price  and  Rainer  Storn,  is  a 
minimization  technique  for  real-valued  problems  which  is  based  on  genetic 
algorithms.  This  project  will  compare  the  performance  of  this  method,  with  the 
performance  of  the  Nelder-Mead  simplex  method.  Both  methods  will  be 
implemented  in  C++  and  called  from  an  existing  MATLAB  implementation. 
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Possible  use  of  Differential  Evolution  in  other  areas 
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Which  one  produces  fastest  solution? 
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Identical  Platfornis  -  Pentium  PC 
Determine  best  method  based  on  speed  of 
detection. 
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Finish  coding  C++  Nelder-Mead  Implementation 
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Finish  coding  C++  DE  Implementation 
Test  and  Compare  Implementations 
Final  Project  Report 
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Efficient  Method  of  Passive  Emitter  Location”  by  Klaus 
Becker  in  the  October  1992  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems. 
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line  segments  and  polygonal  faces. 

Initialization:  Start  with  N  +  1  initial  points  which  give 
an  initial  simplex. 
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■  Stop  when  the  decrease  in  function  value  during  a  step  is  less 
than  some  tolerance 

Restart  at  minimum  . 
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Chapter  10  of  Numerical  Recipes  in  C  (2nd  Edition)  by 
William  H.  Press,  Saul  A.  Teukolsky,  William  T. 
Vetterling  and  Brian  P.  Flannery. 
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Differential  Evolution  is  an  Evolutionary  Strategy 
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them  into  a  special  alphabet. 

Usually  use  a  deterministic  selection  scheme  during  crossover. 

Use  addition  instead  of  symbol-swapping  during  mutation. 
This  makes  incremental  search  easier. 
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member.  Store  these  values  in  a  array. 

V\r 

•  Each  member  of  the  current  population  takes  its  turn  as  the 
target  vector  for  mutation,  crossover,  and  selection 
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Determine  the  difference  between  two  randomly  chosen 
population  members.  xa  -  xb 
Multiply  by  the  scaling  factor  F.  F*(  xa  -  xb) 


w  o 

So  - 
a  s  w 

35k: 

Z  J 
O  ◄ 

P4  > 

H  W 


§  Z 

5  1 

aSs 

sSg 

ix  2  -w 

Cq  g  U 
^  ® 
s  3  ^ 

6  i,  cs 

z  ^  g 

|§| 

gsa 

5  1 1 

t.  Z  S3 

O  bq  >3 

gig 

§gg 

I  g* 

8  S 


8  ^ 

>  -° 

« ? 
o>  1 

C/2  C3 

o  x 

43  w 
O  * 

>,P4 

.  .  'g  + 

^  O 
£  d>  „ 

o  S  11 

CJ  o 


+3  -d 
cd  ^ 


-S  o 
"3  *> 
<  > 


O  *d 

C/3  O 
C/3  *-i 

o  Ph 


4K 

cc  A  fll 
ju  ||  M 


jiH. 

x^/aaotf^ 


!f 


Conduct  D- 1  experiments  starting  at  a  randomly  selected 
parameter. 
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For  each  parameter,  compare  CR  to  a  uniformly  distributed 
random  number  from  0  to  1 . 

If  the  number  >  CR  then  the  trial  vector’s  parameter  comes 
from  the  target. 
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The  cost  of  each  trial  vector  is  compared  to  the  cost  oi  the 
target. 

The  vector  with  the  lowest  cost  is  moved  to  the  next 
generation  array.  If  the  trial  vector  wins,  its  cost  replaces  its 
parent’s  cost  in  the  cost  array. 


B 

*2  g  oo 
£  o2 

d  *  d 

s  S  § 

Is! 


05 

H 

Pm 

W 

U  « 

£  5 

o  3 


™  MH 

eg  5  u 

•<  O  P- 

^  c  ◄ 

^  2  e£ 
U  ^  w 

E2  P 
Z  hJ 

©  ^ 

> 

H  W 


Jgg4 


^  z 

ii* 

III 

S^S 

k5  a.  W 

s  3  ^ 

S  9  fc 

|§§ 

iSg 

IsigS 

>  jJ  j! 

S  £s  ^ 

S  §  ^ 

^  $  g 

5  K  g 

§  £  * 

<j  ^ 


.3  <L> 

S  % 

2  g 

,<U  -G 

5a  h 

•  fH 

Q  * 


3 

**  ^ 
>>  CD 
ccS  vh 

e 

Cd  £ 

4-» 

a  xj 

<D  -— ! 


'S 

o 

CD 

fcs 

g 

5 

1/3 

<D 

•  • 

G 

<D 

s 

O 

O 

•  T— ! 

.s 

o 

<D 

% 

u 

<D 

5-h 

o 

o 

*£> 

£ 

• 

<d 

£h 

CD 

> 

ccS 

a 

o 

a 

o 

+-» 

cd 

•  H 

i  * 

X 

CD 

ccS 

S3  a5 


tj  o 


VL/  LJ 

a  2  5 

©  r^  O 
t>o  H  u 


r-1  o 


"rials  continue  for  a  pre-set  number  of  iterations  or  until  the 
ost  of  the  best  vector  is  less  than  or  equal  to  a  pre-set  value. 


_  s 

"H  2  oo 

O  t >0  On 
Lu  O  On 

U*  r-H 

BS  *  = 

»w  s 

U  C/3  5< 

sa< 

cd  ^  <N 


CZ) 

H 

Ph 

W 

U  5 

S5  5 

o  S 

ug 

So- 
5  &  « 

to  ^  w 

^  z  o 

^  O  Ph 

^  HH  ^ 

^  H  g 

U  ^  w 
H  EJ 
Z  hJ 

o  «! 

Pi  > 

H  W 


«£  Z 

S  | 

U.  a  ■* 

til  g  U 

^  2 

fcj  ij  ^ 

^  S 

l§6 

IIs 

5 1 1 
§g£ 

s  ^  $ 

|§§ 

3  g 


S^^CSL- 


100 


Next  Generation 
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5.1.8  Presentation  by  Dennis  Ludwig 


The  student  briefing  presented  by  Dennis  Ludwig  at  this  meeting  is  reproduced  on  the  next  1 2 
pages. 
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RESEARCH  PROPOSAL  PRESENTATION 

Dennis  A.  Ludwig 

WR-ALC/LNERB 
ALQ-155  System  Engineer 

Background  and  Experience: 

Education:  BS  EE  from  Louisiana  Tech,  1981 

MSA  from  Georgia  College,  1988 
Pursuing  MSEE  from  Mercer,  45  hrs.  completed 
WR/ALC:  Electronic/ Software  Engineer 

Research  Topic: 

Pose  estimation  using  a  Neural  Network 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
(RAPCEval) 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


PROBLEM  STATEMENT 

★  Design  a  Neural  Network  system  to  resolve  the  Pose 
Estimation  Problem  of  a  target/ object  in  three 
dimensions.  That  is  the  estimation  of  location  and 
orientation  of  a  detected  target  in  a  Pose  volume. 

★  Given  an  image  of  an  object  in  a  background,  and  a 
model  that  we  would  like  to  match  to  the  object, 
what  is  the  best  pose  for  the  model?  (Ghee  and 
Doria) 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


PROPOSED  RESEARCH  OBJECTIVE 

Existing:  The  pose  estimation  problem  has  been  solved  for  a 
two-dimensional  image,  involving  three  degrees  of 
freedom  (two  associated  with  position,  one  with 
orientation). 

Improvement:  To  add  a  third  dimension,  which  adds  three 
parameters  (Z,  Phi,  Psi)  to  the  current  two-dimensional 
Pose  Estimation  Neural  Network  System.  This  will 
correspond  to  six  degrees  of  freedom. 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 
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TECHNICAL  SECTION 

3-D  Image  Representation 

.  Three  degrees  of  freedom  for  translation  in  the 
X,  Y  and  Z  directions 

•  Three  degrees  of  freedom  for  rotation  or  orientation 
about  the  X,  Y,  or  Z  axis  (  <x>  ,  w  ,  ©  )  respectively. 

.  Scale  will  be  incorporated  into  a  distance. 

*  Note  that  3-D  Pose  volume  has  6  parameters  describing  it. 


TECHNICAL  SECTION 
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System  Overview 


Image 

(Data) 


Minimum  Distance 
Vector  Calculation 


Pose  Estimation 
Neural  Network 


Model 

- z - 

/ 

r~  ! 

1 _ _ _ 1 

Pose  Adjustment 

L  J 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 

TECHNICAL  SECTION 

Model 

•  An  imitation  of  an  existing  object 

•  Contains  a  priori  knowledge  of  the  feature 

•Is  encoded  in  traditional  data  structure  (two  or 
three  dimensional  array) 

•Referred  to  as  a  synthetic  image 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
( RAPCEval ) 


TECHNICAL  SECTION:  MODEL  CONT. 

Example  of  a  synthetic  image 

A  synthetic  image  encoding  the  knowledge  (perspective 

data)  of  a  road. 
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Data  structure  is  two 
dimensional  array  of 
pixels. 
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TECHNICAL  SECTION:  MODEL  CONT* 


-1, 1,-1 


1,-V1 

■vw 

-1,-1, 1 


Perspective  Inversion  makes  it  difficult  to  see  which  of  the 
opposite  comers  of  the  cube  is  the  nearer. 
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TECHNICAL  SECTION 

Neural  Network 

★  Architecture  (MLFN)  (MLP)  (Backpropagation) 

★  Type  of  Learning  Procedure  (Supervised) 

★  A  set  of  training  patterns  is  used  to  teach  a  NN. 

.  Input/ Output  vectors 
.  Adjust  weights  based  on  output  error 

★  Each  parameter  variation  will  have  an  output  node 
corresponding  to  it. 

★  (  A  X /A  Y  /  A  Z  >  A  <&>  A'L'  ,A<E>) 
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RF  RECEIVER  AND  PROCESSING 
CONCEPTS  EVALUATION  PROGRAM 
(RAPCEval) 


TECHNICAL  SECTION 

Activation  Function 

★  Sigmoid  or  Logistics  function  often  used 

★  Must  be  differentiable  to  use  Backpropagation 


f(x)  =  (l/l+exp(-bx) )  -  0.5 


f(x)  =  f(x)(l-f(x)) 
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Dennis  A.  Ludwig 
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TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 

A  typical  weight  file  for  Five  input  nodes,  Two  hidden  nodes,  and  Six  output 
nodes: 

0(1,0)  0(1,1)  0(1,2)  0(1,3) 

0(2,0)  0(2,1)  0(2,2)  0(2,3) 

0(3,0)  0(3,1)  0(3,2)  0(3,3) 

0(4,0)  0(4,1)  0(4,2)  0(4,3) 

0(5,0)  0(5,1)  0(5,2)  0(5,3) 

0(6,0)  0(6,1)  0(6,2)  0(6,3) 

H(0,0)  H(0,1)  H(0,2)  H(0,3)  H(0,4)  H(0,5) 

H(1,0)  H(l,l)  H(l,2)  H(l,3)  H(l,4)  H(l,5) 
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TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 

Parameters: 

★  Learning  Rate,  can  be  different  for  each  layer 

★  Momentum,  can  be  different  value  for  each  layer 

★  Activation  Function  Temperature 

★  Bias 

★  Number  of  patterns  in  Epoch. 


RF  RECEIVER  AND  PROCESSING  r™™™8 

CONCEPTS  EVALUATION  PROGRAM  update 

( RAPCEval ) 


TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 


_  1.  GATHER  TEST  &  TRAINING  DATA 

START  )  2  CHOOSE  INITIAL  WEIGHTS 


Designing  a 
Backpropagation 
Network 


TRAIN 

SYSTEM 


ADJUST 


N  Training  \ 
Complete^ 


ADJUST 

ARCHITECTURE 


TEST 

SYSTEM 


Good  ^ 
Results?/ 
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TECHNICAL  SECTION:  NEURAL  NETWORK  CONT. 

Data  Collection  for  Training 

★  Set  model  in  pose  volume  and  measure  sampled 
points  from  focal  point. 

★  Re-orient  the  model  and  measure  sampled  points 
from  focal  point. 

★  Determine  the  new  pose  given  those  A_Parairieiers 
(outputs)  and  the  resultant  minimum  distance 
vectors  (inputs) 
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Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 


Training  and  Testing  the  Network 


★  Data  is  divided  into  training  and  testing  sets 


*  Training  occurs  using  collected  data 


*  Trained  NN  is  then  evaluated 


★  Adjustments  made  if  results  are  not  acceptable 


RF  RECEIVER  AND  PROCESSING  ^ef™' 
ill  1  1  1c)  CONCEPTS  EVALUATION  PROGRAM  updat 

_ (RAPCEval) _ = 

PROPOSED  SCHEDULE 

★  Quarter  1 

•  Develop  expertise  in  data  modeling  and  BRL-CAD 
software  (Didn't  use)  , 

.  Develop  expertise  in  neural  network  design 
package  (Didn't  use,  wrote  my  own) 

•  Develop  initial  solution  strategy  for  proposed 
problem 

★  Quarter  2  (Current  Semester) 

•  Develop  and  test  backpropagation  network 

•  Document  results 


Dennis  A.  Ludwig 
MSEE  Program 
UPDATE 
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I  could  get  the  model  to  move  to  the  image  by 
minimizing  an  angle  between  the  light  source 
ray  tracing  the  minimum  and  maximum 
points  of  the  two  images.  The  model  point  is 
well  defined,  but  the  Z  coordinate  will  be 
missing  from  the  target  image. 


Look  at  rotating  an  enclosed  area  to  get  Z 
information. 
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5.1.9  Presentation  by  Tracy  Tillman 


The  student  briefing  presented  by  Tracy  Tillman  at  this  meeting  is  reproduced  on  the  next  30 
pages. 
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Research  Topic: 

HARDWARE  IMPLEMENTATION  for  ADVANCED  PULSE  PROCESSING 


RECEIVER  AND  PROCESSOR 
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DATA  STORED:  RATIO  &  SUM  OF  INTEGERS 
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candidate  countdown  pairs/PRPs 


RF  RECEIVER  AND  PROCESSING 


V 

OO^iOOOLn^-ONroOCMCONN 
C_  ONO‘r)noC0SONOO'^‘CO'r-0 
rS  coONM^^^OLOlOCO^lOCO 
£  C^oqvin^oOOOOT-rOiDOD 

O  ^  ^  ^  ^ _ q  Q  Q  ^ ^ ^ ^ — 


134 


★  91  EQUIVALENT  FP  ADDITIONS 


0 

'^-(NOO^-CSICDCNJ 
CL|  t-LO^-CMUO^'5-O 
rS  LOO'vTCDCDCN'^-CO 
o2  OCOOOOCQCOLOO 
CD^— (N^ 


135 


*59  EQUIVALENT  FP  ADDITIONS 


Xfl 

in 

> 

< 

Z 

< 

m 

& 

P-: 

u 

pa 

< 


c- 


V 


00 

oo 

no 

CO 

ON 

CO 

CO 

o 

NO 

m 

o 

r- 

m 

VO 

m 

oo 

VO 

■^r 

oo 

o 

co 

o 

ON 

(N 

co 

CN 

CN 

CN 

o 

o 

o 

(N 


CO 

UQ 

a 

tLj 


oo  co 

On  — « 


10  CN  f-  CN  i>- 
—  —  ON  OO  VO. 


=3 

CO 


V)  OO  CO  C"*  "sf  '— '  co 

oo  on  — 1  on  oo  c-- 


CO 

t" 

H 

D 

Ch 

5 


H 

D 

Cl, 

H 

D 

O 


-K  -K 


136 


137 


2.2556 


138 


RF  RECEIVER  AND  PROCESSING 


139 


RATIO 


140 


RATIO 


RF  RECEIVER  AND  PROCESSING 


LOOKUP  TABLE 


h-< 

GO 

< 

z 

< 

z 


<  < 

H  (Z 

<  H 

n  w 

M  GO 

w  p 

cn  D 

P  Ph 

s  S  5 

U  ffl  < 
pq  <  w 

<  H  es 

S  2  g 

<  <  Q 
Pi  Pi  pq 

quo 

oo« 

H  H  3 
cn  cn  N 

M  M 

x  x  • 


H 

z 

w 

w 

p 

PO 

H  D 

U  O 

<  x 

W  CO 
pc  W 

w  0  £5 

P  go  P 
PQ  W  H 

H  Z  |> 

P  W  p, 

p  P  O 

^  PQ 

%  D  < 

S  p  cn 
O  U  m 

->  O  p 

O  p  P 

H  o  < 

Q  P  > 
w  w  . . 
X  PQ  H 

S  S  P 

S  D 

Z  2  H 
HH  ^  p 

•  •  o 


*  * 


142 


o 


m 


c f) 

Z 

HH 

< 

P4 

H 

W 

CD 

i_l 

P 

Ph 

H 

P 

Ph 

2 


CO  N  O 


o 

2 


o 

CD 

O 

o 

CM 

CO 

co 

r— l 

r-t 

r“* 

1 

T— ( 

I 

1 

o 

I 

O 

1 

O 

o 

T—< 

r-i 

CM 

On 

m 

r-< 

• 

rH 

• 

• 

cO  nD 
ON  O 


O 

O 


CM  CM 

i  i  i 
CO  N  CO 
CO  O  CO 


O 

2 


143 


166.7  -  213.3 

126.7  -  246.7 


RF  RECEIVER  AND  PROCESSING 


o 

2 

2 

M 

< 

& 

H 

W 
c n 

1 — \ 
)—> 


H 

2 

ID  tNtNNtNNNN^ 

o 

u 


o 

o 

o 


-K 


-K 


144 


652  EQUIVALENT  FP  ADDITIONS 


V  S  g'§L 
£  w  8  < 

r  2  Ci  rs 


05  X 
c/2  O 
a  X 
U  P« 

2  z 

X  o  ^ 

e-5 

0  W 
Zpy 

<  j  Cm 

«  <  <, 
N  >  s 

>  w 

6d  CO 

u  £ 

Sg 

S§ 


5  H 

%  7 

(X  5 

£  8 


NMNNtMNNtMNNl 


sWLLUfP 


IN 

C\ 

O 

— i 

vO 

o 

oo 

CO 

<N 

(N 

CO 

CO 

CO 

uo 

vq 

o 

vq 

o 

7~H 

7H 

7~H 

rH 

t-h 

rH 

7— ( 

rH 

oi 

m. 


M  - 


145 


2.14  2 
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721  EQUIVALENT  FP  ADDITIONS 


.  PULSE  TRAIN  ANALYSIS  WITH  LARGE  TABLES 


5.2  Steering  Committee  Meeting,  August  1998 
5.2.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 
Minutes  August  12,  1998 


A  meeting  of  the  RAPCEval  Advisory  Committee  was  held  in  the  Mercer  Engineering  Research 
Center  Auditorium  at  1:00  p.m.,  August  12,  1998. 

The  meeting  consisted  primarily  of  an  invited  lecture  by  James  P.  Stephens,  Sr..  Dr.  Stevens 
works  in  the  RF  Sensor  Technology  Division,  Electronic  Warfare  Technology  Branch,  of  the  Air 
Force  Research  Laboratory  in  Dayton,  Ohio.  The  lecture  focused  on  the  topic  of  RF  sensor 
technologies  in  the  area  of  C3  Electronic  Combat. 

The  lecture  had  a  communication  countermeasure  emphasis,  focusing  on  advanced  digital  signal 
processing  for  EW.  It  began  with  an  organizational  background  and  objectives,  C2W  trends,  and 
vision  from  the  perspective  of  SNRW. 

A  brief  tutorial  was  provided  to  explain  the  motivation  of  the  use  of  time-frequency  distributions, 
wavelets,  and  cyclostationary  signal  processing.  Current  research  activities  of  SNRW  in  all  areas 
related  to  C3  Electronic  Combat  were  introduced.  These  include  such  areas  as  featureless 
waveform  countermeasures,  time-frequency  analysis,  frequency-hopping  prediction,  GPS  ECM 
and  ECCM,  and  covert  waveform  detection  and  jamming.  In-house  activities  include  time- 
difference  of  arrival  using  cyclostationarity,  frequency  hopping  synthesizer  development, 
spectral  correlation  analysis,  and  transform  domain  waveform  experimentation. 
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5.2.2  Attendance  Roster 


Name 

(last) 

Name 

(First) 

Organization 

;  Phone 

C'\-  li' !  '$$$$ 

E-mail  | 

Bahl 

Erik 

MERC 

912-953-6800 

ebahl@merc.mercer.edu 

Bass 

Charles 

MERC 

912-953-6800 

cbass@merc.mercer.edu 

Bass 

Tom 

MERC 

912-953-6800 

tbass@merc.mercer.edu 

Black 

Joseph  K. 

MERC 

912-953-6800 

jblack@merc.mercer.edu 

Bryant 

Peter 

MERC 

912-953-6800 

pbryant@merc.mercer.edu 

Campbell 

Mark 

WR-ALC/LYSFR 

912-926-5673 

mcamp- 

bel@avionics2.robins.af.mil 

Carter 

Glenn 

MERC 

912-953-6800 

gcarter@merc.mercer.edu 

DePuy 

George 

MERC 

912-953-6800 

gdepuy@merc.mercer.edu 

Finnigan 

Skip 

MERC 

912-953-6800 

sfmnigan@merc.mercer.edu 

Guest 

Robert 

MERC 

912-953-6800 

rguest@merc.mercer.edu 

Harris 

Rhondia 

MERC 

912-953-6800 

rharris@merc.mercer.edu 

Hillhouse 

Ed 

ITC 

912-922-8106 

hillhouse@itc-2.com 

Jones 

Houston 

TASC 

912-926-5458 

jhjones@hom.net 

Kelley 

Joseph  F. 

WR-ALC/LNIE 

912-926-4482 

kelley.joe@ec.robins.af.mil 

May 

Mark 

MERC 

912-953-6800 

mmay@merc.mercer.edu 

McFarland 

Chris 

MERC 

912-953-6800 

cmcfarland@merc.  mercer, 
edu 

Messer 

Bob 

MERC 

912-953-6800 

bmesser@merc.mercer.edu 

Miller 

David 

WR-ALC/LNEV 

912-926-2994 

miller@ec.robins.af.mil 

Mitchell 

Ray 

MERC 

912-953-6800 

rmitchell@merc.mercer.edu 

Moody 

Douglas 

SAIC 

912-918-2911 

moodyd@saic.com 

Moody 

Douglas 

SAIC 

912-918-2911 

moddyd@saic.com 

O’Dowd 

Bill 

Ball 

912-922-4363 

bodowd@f1 5aisf.robins.af.mil 

Oliver 

R.P 

WR-ALC/LNERT 

912-926-2588 

oliver@ec.robins.af.rn.il 

Pequignot 

Nick 

AFRL/SNRD 

937-255-6127 

pequignot@sensors.af.mil 

Potts 

Marvin 

AFRL/SNRW 

937-255-5579 

pottsmn@sensors.wpafb.af. 

mil 

Smith 

Jermaine 

MERC 

912-953-6800 

jsmith@merc.mercer.edu 

Smith 

Jerome 

WR-ALC/LNEV 

912-926-2994 

smithjm@ec.robins.af.mil 

Thorstens 

Gary 

Ball 

912-922-4363 

gthorste@ball.com 

Tillman 

Tracy 

MERC 

912-953-6800 

ttillman@merc.mercer.edu 

Wimpy 

Dale 

WR-ALC/LYSFR 

912-926-9529 

dwimpy@f1 5aisf.robins.af.mil 
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5.2.3  Overview  of  the  Program  (Dr.  Bass) 


The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  1 1  pages. 
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Committee  to  be  uselul  to  tne  Air  t  orce  ana  is 
certified  by  the  University  as  having  academic 
merit 
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Communication  of  Time  and  Time  Measurements”,  unclassified, 
available  from  MERC 
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★  Kirk  Wright,  “Object  Oriented  Modeling 
172”,  classified,  draft  in  edit  at  AFRL 


£ 

O 


C/5 


w 

P4 

Ph 


C/5 

* 

< 

O 

O 

H 


•\ 


£ 


•b 


158 


159 


experimentation.  This  90-minute  briefing  is  intended  as  a  technical  overview  of  the 
current  activities  of  SNRW  and  should  be  of  interest  to  technical  managers, 
researchers,  faculty,  and  engineering  students. 
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in  a  three-state  area  with  prior  assignments  in  Detroit  and  Norfolk.  Since  Jan  1982,  Mr.  Stephens  has  been 
employed  part-time  as  an  instructor  in  the  Electrical  Engineering  Technology  curriculum,  Associate  Degree 
program,  at  Miami  University  where  he  teaches  a  variety  of  courses  in  electronic  and  microprocessor 
technology.  He  has  authored  or  co-authored  more  than  15  publications. 


5.2.4  Presentation  by  Jim  Stephens,  AFRL/SNRW,  Dayton 


A  guest  lecture  presented  by  James  P.  Stephens,  Sr.  at  this  meeting  was  described  by  a  synopsis 
in  the  overview  of  the  previous  section.  Copies  of  this  presentation  are  available  on  request  from 
the  RAPCEval  Program  Manager  at  Mercer  Engineering  Research  Center. 
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5.3  Steering  Committee,  November  1998 
5.3.1  Meeting  Minutes 


EW  Receiver  and  Processing  Concepts  Evaluation  Program 
Program  Research  Standards  Committee  Meeting 

Minutes  9  November,  1998  _ 


The  RAPCEval  Steering  Committee  held  a  meeting  at  MERC  on  9  November,  1998  at  1 :30 
p.m.  Dr.  David  Barwick,  MERC  Executive  Director,  kicked  off  the  meeting  with  a  welcome. 

Dr.  Tom  Bass,  MERC  RAPCEval  Program  Director,  then  presented  a  “State  of 
RAPCEval”  slide  overview  presentation. 

Dr.  Bass  then  introduced  the  three  students  who  were  to  give  presentations  on  their  research 
topics.  Those  students  and  their  topics  were: 

♦  Mr  Houston  Jones  -  “Evaluation  of  Reed-Solomon  Codes  for  CDMA  Systems" 

♦  Mr  Mark  Napier  -  “Civil  IFF  Reed-Solomon  Code  Application” 

♦  Mr.  Randy  Ford  -  “Passive  Location  via  Evolutionary  Genetic  Algorithms” 

The  next  agenda  item  was  a  period  of  discussion  and  suggestions  about  the  projects.  Ches 
Rehberg  (during  discussions  and  via  a  later  e-mail)  pointed  out  the  need  for  Air  Force  personnel 
adept  in  the  topics  being  presented  to  attend  the  meetings.  For  example,  two  of  the  topics  were 
on  communications  (robustivity  via  Reed-Solomon  code  implementations).  Ches  suggested  such 
participants  for  future  meetings. 

Phil  Oliver  stated  that  data  links  are  a  big  topic  in  the  F-15  community.  He  suggested  that 
these  communication  robustivity  research  topics  may  have  an  application  in  that  area. 

Several  attendees  remarked  that  the  student  presentations  appeared  higher  in  quality  on  this 
occasion.  Some  questions  were  raised  about  the  applicability  of  the  communication  topics  to  the 
original  RAPCEval  charge. 

The  meeting  was  then  adjourned. 
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5.3.2  Agenda 


Agenda 

RAPCEval 

program 


STEERING  COMMITTEE 
MEETING 

98  November  9  —  1:30  PM  to  4:15  PM 
Conference  Room,  Mercer  Engineering  Research  Center 
135  Osigian  Boulevard,  Warner  Robins,  GA 


Meeting 
called  by: 

Nicholas  Pequignot,  AFRL/SNRP, 

Air  Force  Research  Laboratory 
Program  Manager 

Facilitator:  Dr.  Tom  Bass 

Committee 

Members: 


Email 

addresses: 


AF  Research  Laboratory 

Mr  Nicholas  A  Pequignot 
Mr  Emil  R  Martinsek 
Mr  Norman  A  Toto 
Dr  Duane  A  Warner 
Mr  Paul  J  Westcott 

pequigna@sensors.wpafb.af.mil 

martinser@sensors.wpafb.af.mil 

totona@sensors.wpafb.af.mil 

wamerda@sensors.wpafb.af.mil 

westcopj@sensors.wpafb.af.mil 


Robins  AFB 


Mr  Steve  Strawn 
Mr  John  LaVecchia 
Mr  Phil  Oliver 
Mr  Ches  Rehberg 
Mr.  Larry  Sheets 

lnerr@ewir-wr.robins.af.mil 

lavecchia@ewir- 

wr.robins.af.mil 

oliver@ewir-wr.robins.af.mil 

rehberg@ewir-wr.robins.af.mil 

sheets@ewir-wr.robins.af.mil 


Mercer  University 
Dr  David  Barwick 
Dr  Tom  Bass 
Dr  Aaron  Collins 
Dr  Benham  Kamali 
Dr  Paul  MacNeil 

dbarwick@merc.mercer.edu 
tbass@merc .  mercer .  edu 
coIIins_as@mercer.edu 
macneil_pe@mercer.peachnet.edu 
kamali_b@mercer.peachnet.edu 
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Meeting  Agenda,  continued,,, 


Schedule 

Greetings 

Dr  David  Barwick 

1:30  PM -1:35  PM 

Meeting  Overview 

Dr  Tom  Bass 

1:35  PM -1:45  PM 

Student  Proposals 

Students: 

Evaluation  of  Reed- 
Solomon  Codes  for 
CDMA  Systems 

Mr  Houston  Jones 

1:45  PM -2: 15:  PM 

Civil  IFF  Reed- 
Solomon  Code  Ap¬ 
plication 

Mr  Mark  Napier 

2:15  PM -2:45  PM 

Passive  Location  via 
Evolutionary  Ge¬ 
netic  Algorithms 

Mr.  Randy  Ford 

2:45  PM- 3:15:  PM 

Discussions  &  New 
Business 

Dr.  Tom  Bass 

3:15  PM -3:45  PM 

Adjourn 

4:15  PM 
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5.3.3  Attendance  Roster 

The  attendees  at  this  meeting  are  listed  here: 


Name 

Organization 

Phone 

Email  Address 

1 

Tom  Bass 

MERC 

912-953-6800 

bass_wt@mercer.edu 

2 

Dave  Barwick 

MERC 

912-953-6800 

barwick_dt@mercer.edu 

3 

Ches  Rehberg 

WR-ALC 

/LNEX 

912-936-4525 

rehberg@ec  .robins  .af.mil 

4 

Charles  Bass 

MERC 

912-953-6800 

bass_cd@mercer.edu 

5 

RP  (Phil) 

Oliver 

WR-ALC 

/LNERT 

912-926-2588 

oliver@ec.robins.af.mil 

6 

Steve  Strawn 

WR-ALC 

/LNERR 

912-926-6435 

strawn@ec.robins.af.mil 

7 

Tony  White 

AFRL 

/SWRP 

937-255-6127, 
Ext.  4236 

whiteaw@sensors.wpafb.af.mil 

8 

Nicholas 

Pequignot 

AFRL 

/SNRP 

937-255-6127, 
Ext.  4235 

pequigna@sensors.wpafb.af.mil 

9 

Mark  Napier 

Scientific  Atlanta 

770-903-6980 

mark.napier@subasic.sciatl.com 

10 

Randy  Ford 

WR-ALC 

/LYSBC 

912-926-0423 

fordra@avionics.robins.af  mil 

11 

Houston  Jones 

WR-ALC 

/QLYM 

912-926-5458 

jhj  ones@hom  .net 

12 

Behnam 

Kamali 

Mercer 

912-752-2415 

kamali_b@mercer.edu 

13 

Paul  McNeil 

Mercer 

912-752-2185 

macneil_pe@mercer.edu 

14 

Joseph  K. 

Black 

MERC 

912-953-6800 

jblack@merc.mercer.edu 

15 

Ray  Mitchell 

MERC 

912-953-6800 

rmitchell@merc.mercer.edu 

16 

Skip  Finnigan 

MERC 

912-953-6800 

sfinnigan@merc.mercer.edu 
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5.3.4  Overview  of  the  Program  (Dr.  Bass) 


The  Overview  Briefing  of  the  RAPCEval  Program  as  presented  at  this  meeting  is  reproduced  on 
the  next  7  pages. 


166 


167 


168 


AH  research  has  been  determined  by  the  committee 
to  be  useful  to  the  Air  Force,  as  well  as  having  merit 
for  graduate  research 
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LaVecchia,  Mr.  Phil  Oliver,  Mr.  Ches  Rehburg,  Mr. 
Larry  Sheets 

*  MERCER  UNIVERSITY:  Dr.  Aaron  Collins,  Dr. 
Benham  Kamali,  Dr.  Paul  MacNeil 
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Solomon  Codes”  to  be  included  in  forthcoming  unclassified 
report 

Mark  Campbell,  “Auto-Regressive  Spectral  Analysis  - 
EW  Applications”,  unclassified  report  available  from  MERC 
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★  Joseph  Kelley,  “MultiGroup  Simultaneous  RAD 
Parameter  Selection”,  WL-TR-97-1094  (classified). 
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EW  RECEIVER  AND  PROCESSING 
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5.3.5  Presentation  By  Mark  Napier 

The  student  briefing  presented  by  Mark  Napier  at  this  meeting  is  reproduced  on  the  next  1 5 
pages. 
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Research  Topic: 

APPLICATION  OF  REED-SOLOMON  ENCODING  TO  IMPROVE  PROPOSED 
COLLISION  AVOIDANCE  SYSTEM  BASED  ON  CIVILIAN  ATCRBS 
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reliability  improvements  to  be  realized  from  the  RS  code.  A  decoder  will  be 
designed,  implemented  and  tested  in  Verilog. 


•  Similar  encoding  could  be  used  in  Military  IFF 

•  Military  aircraft  could  participate  in  the  civilian 
system  without  revealing  aircraft  ID 
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Design  the  Decoder  in  Verilog 
Verify  with  Verilog  Testbench 
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Model  without  Erasure  Information 
Model  With  Erasure  Information 


.2  d  oo 

ON 

2  P  ^ 

(C  r/  h 

S  S  - 


O  <3 

5T  P4 

r  * 


w  w 
u  ^ 


O  2 
P4  O 


Cd  H  ^ 

o^« 

z3£ 

<j  j  3 
P4  ^  ^ 

w  >  C 
>  w 

E  52 
u  fc 

w  S3 

tt  S 

fe  z 

K  O 

U 


I 


r® 


i  1° 


Q 

w 

CD 

O  LO 

OQ 

O  o' 
oc  o 

Oh  b- 
tq  ^ 

9  3 

Bj  ^ 

gfe  w 

ou  hJ 

8  2  D 

o°D 

§Sw 

g3£ 
8  *  u 

25  t/l 
w  w  _ 

2;  Q 
o£  W 

5  ^  CD 

O  tq  O 

«  ^  O 
8§* 
3§  ^ 

O  £ 

O  cd 


E-h  -J 
rr1  1-q 

0  o 

Cq 

Oh 


ci 


*-<  <D 

T3 


CD  >rH 

os  ,2 
^  Pin 
<3  Ph 


0) 

i _ \ 

03 

<L> 

?H 

__f 

CD 

03 

CD 

•  rH 

CD 

-M 

Jh 

OS 

pH 

03 

X 

03 

Vh 

<D 

CD 

C 

*4-> 

•  rH 

60 

03 

i  1 

C 

<D 

X 

•rH 

u 

?H 

03 

CD 

CD 

c 

60 
•  fH 

cd 

•pH 

> 

O 

u 

P-. 

<L> 


60 

.a 

T3 

03 

M-l 


G  . 
fH  60 
CD  £ 
£j  «Jh 

5  T$ 


£  O 

^  CD 


.5  $ 

CD  P4 


Develop  and  test  decoder  for  RS(31,23)  code. 
Document  results  and  submit  project  report. 


s 

O  <1 

a  o 

CD  X 

cn  y 

W  K 

(J  pH 

O  2  _ 

g2  -a 

<  j  3 

X  <  ^ 

w  >  C 
>  w 


Q 

tq 

co 

O 

Bh 

O 


to 

CQ 


tt$  u 

&H  t-H 

tq  ^ 

O  ^ 
c<  ^ 
Bh  tj 


O 

E-. 

O 

£ 


u 

£ 

O 

Q 

tq 


Q  2 
O  ^ 

o  « 

£  s 

tq  tq 
E— i 
co 

>H 

CO 


£ 

o 


O 

hJ 

O 

CO 


tq 

Q  2 
tq  ^ 
tq  O 

C<  ^ 


tq 

o 

£ 

g 

s 

•~H 

►q 

tq 

Cq 


£ 

o 

I— H 

CO 

I-H 

-q 

►q 

O 

U 


Z 

o 

HH 

H 

< 

H 

Z 

W 

CD 

W 

C* 

Ph 

hJ 

< 

u 


Z 

EC 

u 

w 

H 


j-i 

VH  43 

°  'g 

CQ 

•g  'o 
£  is 

d  g 
o  o 

*3  U 
^  u 


a  • 

t_i 


.2 

a 


cQ 

5h 

H 


pq  j=j 

&  i 

HH  tO 

S  PQ 
g  C* 

^  y 

«*  H 
a 

O  pQ 

’-O  ^ 
©  2 
CD  ^ 

Ctf  X) 

©  C 
b  a 

C3  q 

2  bb 
x  o 

c  Lh 
Oh'S 

3  | 

03  a 

Vq  © 
q->  h—> 

a 

C^,  -i-H 

a  22 

o 


bo 
d 

*  ?H 

CD 

a 

4-> 

o 

C  a 
<£  o 

H-q  ,^h 
03  CD 

i-H  *f-4 

o  ;z5 
.b  2 
*3  o 

o  G 

q-q  h— > 

•  fH 

£  £ 
*5  00 

eg  vh 

3 -8 

£  § 

M  ° 
•*H  J_< 

*  fH 

03 

i-q 

O 

‘S? 


a 

^o 

o3 


i-H 

*2 

c 


T3 

3 

o 

1  f 

Cl  <D 
Jh  h 

a  c« 
t>  >> 

T3  CO 

2  bo 

a  -s 

I  § 

►>*>  ^ 
Vh 

<u  JO 

>  o 

fq 

cQ  ra 

^  5 

cQ  O 

-s  o 
a  to 

<U  Dh 

tS  O 


CD 

© 


CD 

© 

> 


o  > 

£  2 
03  Bh 

fc  »2 

C  < 

<d  n 

bo  H 


bS) 

c 

-i-H 

CD 

a 


S 

© 


ctf 

a 

cd 


o 


X5 

a 

03 

© 

CD 

a 


(D 


a 

3 

CD 

>> 

CD 


CD  © 
eCS  O 
H  &H 


CD 

© 

O 

•H 

c: 

<D 

&o 

'o 

H 


_  d 

CO  O 

>-,  o 
CO  o 

g  £ 

8  2 
c^S  w 

CD  i-H 
Jx  *1-H 

CQ  o3 


(/3 

d 

<D 

CL  _ 
X  W3 
<D  >% 

d  “ 

CQ  O 

.22  e 

a  * 

g  3 
x>  o 


d 

o 

•  fH 

c^S 


> 
03 

a 
o 

S  2 

i-q 

o  ^a 
o 

a  o 


"d 
© 
q-q 

3 

_  £ 
»-q  i-q 


© 

O 

a 

o3 

tj 


o 

> 

03 


CD 


rfD 

3 

c3 

■  f-h 

"© 

i-H 

2 

W) 

•  fH 

fC 

X5 

a 

03 

© 

> 

•  f-H 

CD 

C 

© 

Pin 

X 

© 

a 

•  t-q 

<u 


181 


<D 

c/3 


O  _ 

o  Q 

I  hh 

G  * 

o  > 


o  a 

P 

S  w 
^  0) 
£  33 
O  § 
33  C 

a  £ 

o  2 

P  g 

«2  -2 
G  H 
G  G 

h  > 


mh  3 

ffi  .ts 

5  43 

Hn  Cd 

®  44 

8  | 

^  g 

g  cr 

O  C/3 

'i  u 

o 

S 

C/3  i-H 

CD 

4"i  i)  /''""‘N 

G  <U 
50  33 

2  ° 

Vh  ^ 

£  44 
G  £ 


c/5  — 

c«  _i< 
*-  ^ 
G  a> 

2  33 

fc:  o 

3  g 

o  1-1 

<U  ’S 

43 

H  £ 


.  33 
>,  £ 

5  •- 

2  s 

*G  2 

P-  5 

n  2 

a  w 

C/3  O 
**  c"* 

<3  ^ 
-o  & 

G  JT 
O  ^3 

P 

2  o 


50 

2  c 


*8  % 
>  3 
•^h  cl 
Q  *7 
o  ^ 
<D  u 

lL  c« 

« = 
s - 

<u 

3  c/5 

P  ■*— “ 

CT  *-J 
<D  <3 

a  ^ 
o  o 
•a  >o 

Cd  (N 
50  cm 

2  ° 
c  t; 
<d  ? 


<D  O 

■S  s 

G  £ 

©  O 
50  P 

•S  M 

5  <u 

§  G. 

S’  2 

-a  a 

2  M 
O  4-> 

S?  S 

C/3  c/3 

c/3  r* 

2  § 

6  is 


N  | 

ffi  PQ 

Is 

s° 

c/5  P 

G  O 
•2  >> 
8  ;| 

°  •  3 

c/3  <-h 

a  g 

D.  “ 
m  33 

2  O 
n>  43 
oo 

£  p 

g£ 

i—L  w 


73  ^ 
^  g  0) 

o  o 
o  ^ 

S  ^  12 

c3  re 

*-■  ±3  g 

g  2  ° 

4-4  — ^  Lh 

•h  C  o) 

£  o  > 

W5  *-m> 

M  *M  <U 

S  73  O 

2  'O  ox 

i3  es  *^ 

o  G  -2 

— <  £) 
>>  D  *  rt 

WX-O  cd 

o  3  Oh 

O  -3  S 

G  g*  O 
x:  g  0 
O  o  w 
0)  -Sh  ^ 

*-  h  js 

h  <u  2s 

-§  g  g 

G  O  w 
O  5S  cm i 
p  ^  cd 

G  X5  2 

2  >  ^ 
G  50  G 
a)  r*  cd 

y  -2  73 

g  G  <0 

^  ^  on 

<D  ^  G 

C/5  **H 

3  O  >> 
T3  ^-h  W 

g  >  % 

^  ^  HH 

, — i  cd  M 

ZL  G  o5 

D  o  a) 

E2  S 
«  x  o 

43  O  p 

O  P  ^ 

“  CO  2 
33  n, 

8  3  ^ 

< 

O  G  ^ 

>  <D 

5“  fe  •§ 

^  c  c 

•  i—i  c 


C/3 

>.  50 

C/5  O 


p  G 
43  •’“• 
H  ^ 

^  o 

-S 

2  ,5P 


<a  m> 
43  G 

Vh 

G  G 


2  G 


c3  ’■^ 

w  OJO 

Vh  •*-< 

a)  — 

«  3 

50  <D 
G  43 

•  M  -M 

^  o 

O  ^ 
43  co 
co  G 

>>  ° 

s  g1 

33 

c/5  G 

43  c 
33  cd 

4-4  *n 

•  i-H  rv 

P  P 
4^  43 
O  50 
O  -d 
o  hJ 


G 

P 


G 


Position  Modulation)  scheme  with  a  1  Mbit/s  rate.  A  “1”  is  defined  to  be  a  0.5us  burst 
followed  by  0.5us  of  off  time.  A  “0”  is  defined  to  be  0.5us  of  off  time  followed  by  a  0.5us 
burst.  The  message  is  proceeded  by  a  8us  sync  pulse.  Either  56  (single  length)  or  112  (double 
length)  bits  of  data  follow. 


H  oo 

«  co  2 

-9^ 

fg  f/  H 
^  pH  T_H 


o  g 

^  P5 

s  o 

cn  g 

</>  5* 

W  K 
U  Ph 

02^ 

go-5 

CL(  HH  W 

Q  ^  W 

<3  5  & 
^  ^ 
>  w 

S  52 

□  fa 

w  fig 

*  u 

b  z 

*  o 

u 


|5|  ]|  '""I  i  pj 


U  .S' 

pj  X 

E  ^  « 

W  oj  M 
—  Qh  03 

O  o  +-> 

gs2 

»a» 

fc  s  s 

o  -*-> 

^  2  x 

fa  O  X 

2  ^  ^ 
Vh  fH  rrf 

tq  a  ™ 

-O  8  « 

fa  X  2 


S-H 

o  <D 

r  i. 


-i  *-*  L 

®3  © 

§  “  ■§ 

^  *  s 

(D 

>  co  m 

52  2  ^ 

-^3  X  ^ 

+2  <l>  cn 

•  fh  Co  s— ''/ 

X  C  00 
OwP^ 

S-H 

<D  W)  ° 

M  C  t 
a  X  II 

^  X  w 
C/3  O  4-> 
<D  O 
rj  C  £ 
C  (D  ^ 

■5  «  C3 
£?  |  x 


<D  ’o 

3  °? 
J’S 

T3  D 
<D  f* 

■5  w> 

Lj  $-4 

5?  •*-( 

O  C/5 

ti  3 


D  .0 
PD  co 

C  ^ 

S  2 

I  & 

^  X 

£  • - 
r?  £ 


o  a 
E  "*  V, 

fr _ r  <L> 

S  'o  w> 

h  1)  C3 

M  a  § 
<3  S3  In 

S  -8  5 


<u 

fa 

O  cS 


cd 

O  ts 

r-H 

4_> 

fa  X 

r-H 

o 

O 

V# 

V# 

2 

Vh 

a> 

«  o 

73 

C3 

o 

o 

ts  ^ 

5-h  ^ 

•> 

c 

s  . 

0 

cd 

o 

*s  s 

0 

VJ 

(D 

<u 

73 

<3  x 

0 

O 

>  X 

C3 

O 

C/5 

'53  £ 

O  J2 

2  §* 

•  fH 
00 

•  fH 

fC 

r 

H 

u  u 

3 

. 

e 

cd 

<u 

^  0 

73 

>,  xi 

cd 

o 

fO  0 

> 

O 

^  (U 
73 

cd 

<D 

co 

X  £ 


|  § 

£  % 
Cd  & 

&  *$ 

s  * 

a  a 


<D  ^ 

73  o 
>  2 
2  a 

&  C 

<D  « 
X>  4> 

a  5 
S  g> 

C 

.2  3 

s  § 
£  -§ 
o  £> 

S  > 

-fa 

2  o 
2  £ 


73 

Vh 

2 

(U 

O 

(D 

S3 

•  fH 

fcs 

cd 

<D 

•  fH 

4-^ 

fO 

-t— > 

00 

O 

0 

co 

C/5 

•i-H 

35 

fO 

< 

Cfl  Y  - 
X  _, 

o  -o  'H 

.2  3  g 

e  -s  £ 

fa  co  « 
«2  «  ^ 
•S  c  « 
2  §  o 
o  cr  ft 

M  <D  -a 
ccJ  &o 

£  u  2 

co  -O 
W)  1)  ^ 
C  *C  Vh 
.0+^0 
e  r-  4h 

*3  ’S  co 
5  *3  TJ 

O  _,  O 
<L>  ^ 

4-h  .  Sh 

S<U  X) 

_  o  ^ 
«  ^  2 
^3  .2 


co  <L> 
O  % 
°  2 

co  a 
r±4  <D 

O  K. 

Gt 


O  £ 

3  W» 
^  3 


3  co 
-H  O 
g  Oh 
*-i  O 

2  ° 
«  Oh 


!-i  ft 

2  ° 
C3  Oh 

<D  7/ 

b£)  2 

M  S 

,0  D 
4h  33 
•H  O 
UZ  03 

s « 

a)  W 
X  tin 
3  <D 
D 

x  H 

T3 

32  « 

5  ° 
o  e 

5  ^ 

E  o 

<D  > 
to  ^ 
>-»  c3 
“  O 

V,  '33 

O  -t-h 
C/5  ^ 

Oh  8 

2  g 

Oh  £ 


c  o 

.2  3 

C/5  C/5 

5  <U 

o  > 

s  ^ 

o  o 

O  .0) 
C3 

kh  <U 


Vh 

(D 

T3 

O 

O 

(D 

^  d 


C/5 

0^  g 

D  X 
*  ^ 


c^  m 
^  in 
1-4 

i 


32  >>  X 


‘r2  ’-' 

-2  o 
.2  x 

13  M 

0  co 
3h  <U 

s  § 

<L>  p 

S«  cd 
>>  rj 
00  “ 


2  S 

3?  C/5 

> 

O  03 


s 

o  ^ 

y  & 

K  O 

CD  Q 
CD  y 

W  p 

(J  Ph 

02^ 

es-a 

5  =  2 

^  H  < 
Pi  5< 
w  >■  r2 
>  tu 

S  S£ 

U  fc 

W  m 

« s 

ft  z 

K  O 

u 


a  £ 

2  5 
O  X 
O  O 

<?  00~ 
X)  On 

0)  r  >> 

<u  u 
£*  H 
> 
°  PJ 

c/3  rl 

<D  W 

•z3  pq 


m-h 

•§  ® 
Q,  V3 

U  .2 

c  ^ 

5  <u 
.2  o 
+■>  o 
O  O 

<D  ?-i 

fc! 

o  „  - 
O  m 

5-1  5 
o  o 


Ui 

W 


a 

_o 

"2 

3 

2 

2 

o 

CJ 

I 

JL> 

3  r-' 
O  5 

s <? 

•2  ^ 
on  cn 


C 

o 


&. 

CX 


3  „ 

o  oo 
.3  On 
CX  On 
CX  ’-' 

c3 


< 

J3 


GO 

0) 

TJ 

O 


c3 

T3 

?3 

C 

c3 


u  u 


0) 

o 


§ 

& 


r\ 

C/3 

e 

o 

•  — H 
•*-» 
c3 

’a, 

ex 

< 


o 

s 

C/3 


o 

o 


Q 

^  • 
.2  oo 

_)  On 

M  r— < 


od  — r 

r— '  '«$ 

£2,  ffi 


184 


Mode  C  Packet  -  Identical  to  Mode  A  packet.  Altitude  encoded  on  10  bits  of  the  digit  values. 
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Schedule  graduation  (campus) 
Project  committee  defense 
Submit  project  report  to  university 


5.3.6  Presentation  by  Houston  Jones 


The  student  briefing  presented  by  Houston  Jones  at  this  meeting  is  reproduced  on  the  next  15 
pages. 
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SOLOMON  CODES  IN  WIRELESS  CODE  DIVISION  MULTIPLE  ACCESS 
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Applicability  to  the  USAF 

|  •  Digital  cellular  communications  are  important  and  valuable  to 
military  as  well  as  civilian  users .  Increases  the  knowledge  base 
in  digital  mobile  communications 
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Determine  the  feasibility  of  using  Reed-So 
(RS)  codes  to  reduce  system  complexity  in 
CDMA  mobile  communications  systems 
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Optimistic:  Improved  performance  and  reduced 
complexity 
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Study  to  determine  the  feasibility  of  using  RS 
code  families  for  single  burst  error  correction  in 
Rayleigh  fading  channel 
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Error  correction  capability 
Frame  error  rate;  bit  error  rate 
Delay 

System  complexity 
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In  1993  a  cellular  system  based  on  code  division  multiple  access 
(CDMA)  was  developed.  This  system  uses  direct  sequence  spread 
spectrum  techniques  requiring  less  power  and  allowing  for  more 
users-than  the  other  designs: . . . — . . . - . 
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Forward  CDMA  channel  modulation  process. 
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Submit  thesis  to  university 


5.3.7  Presentation  by  Randy  Ford 


The  student  briefing  presented  by  Randy  Ford  at  this  meeting  is  reproduced  on  the  next  47  pages. 
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Research  Topic: 

COMPARISION  OF  DIFFERENTIAL  EVOLUTION  TO  THE  SIMPLEX  METHOD  IN 
OPTIMIZATION  DURING  PASSIVE  EMITTER  LOCATION 
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Differential  Evolution,  developed  by  Kenneth  Price  and  Rainer  Stom,  is  a 
minimization  technique  for  real- valued  problems  which  is  based  on  genetic  algorithms. 
This  project  will  compare  the  performance  of  this  method,  with  the  performance  of  the 
Nelder-Mead  simplex  method.  Both  methods  will  be  implemented  in  C++  and  called 
from  an  existing  MATLAB  implementation. 


I 

I 

1 

S 

8 

I 


I 

5 

o 

I 


s 

to 

3 


s  £ 


I 

I 

I 

I 

g 


s 

I 

N 


K 

% 

£ 


S 


1/1 


m 


m 

<D 


m\ 


209 


210 


Identical  Platforms  -  Pentium  PC 

Determine  best  method  based  on  speed  of  detection. 
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★  Test  both  implementations  (cont) 

•  Compare  accuracy  of  solutions  obtained. 
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x,y  plane. 
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Becker  in  the  October  1992  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems. 
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an  initial  simplex. 
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Low  |  All  Dimensions 
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Algorithm  that  is  especially  useful  in  minimizing 
functions  with  real  variables  and  many  local  minima. 

Differential  Evolution  is  an  Evolutionary  Strategy 
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Each  member  of  the  current  population  takes  its  turn  as  the 
target  vector  for  mutation,  crossover,  and  selection 
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population  members.  xa  -  xb 
Multiply  by  the  scaling  factor  F.  F*(  xa  -  xb) 
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parameter. 
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The  last  parameter  always  comes  from  xc'  in  order  to  make  sure 
that  the  trial  vector  differs  from  the  target. 
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Trials  continue  for  a  pre-set  number  of  iterations  or  until  the 
cost  of  the  best  vector  is  less  than  or  equal  to  a  pre-set  value. 
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Next  Generation 
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Calls  external  C++  dlls  called  MEX  files  for 
minimization 

Two  versions  -  nmloc.m  and  deloc.m 
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Both  read  in  data  from  same  external  file  data.m 
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data.m  floatf.exe  nmc275.dll 


% 

o 


m 


O 


O 

W 

H 


•  •  f*  V 

-  t!  S  « 

C  S3  O  m 

"  a«  fc  B 


S  ?  «  ®  20, 

«  5  g  fc  .2  W  a 

(/)  O  2  0>  'O  AS 


O* 

0> 

> 


236 


data.m  floatf.exe  dec95.dll 
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Searches  the  error  space  by  calling  amoeba. 
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Searches  the  error  space  by  calling  de_minimizer. 
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and  more  accurate  in  only  5  of  the  12  test  sets. 
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Further  study  of  alternate  crossover  techniques  may 
make  DE  faster  and  more  accurate 
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5.4  Thesis  Report  by  Ron  Brinkley 

A  version  with  standard  report  style  has  been  generated  from  the  Master’s  thesis  of  Ron 
Brinkley.  This  is  reproduced  here  on  the  next  83  pages. 
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1  Summary 


Communication  and  storage  systems,  whether  they  are  magnetic  or  optical  recording  de¬ 
vices,  satellite  or  cellular  mobile  communication,  or  even  military  communication  sys¬ 
tems  in  hostile  jamming  environments,  are  faced  with  a  need  to  correct  burst  errors.  It  is 
the  intent  of  this  research  work  to  explore  the  efficiency  of  powerful  burst  error- 
correcting  Reed-Solomon  (RS)  codes.  First,  an  explanation  as  to  the  cause  of  burst  er¬ 
rors,  and  a  definition  of  burst  errors  is  presented.  This  explanation  is  followed  by  an  in¬ 
troduction  to  the  finite  field  mathematics  necessary  in  understanding  RS  codes.  An  in¬ 
troduction  to  RS  codes  and  their  ability  to  combat  burst  errors  is  given.  Subsequently  an 
overview  of  the  most  widely  used  interleaving  techniques  is  presented  to  complete  the 
preliminary  portion  of  this  thesis. 

The  guaranteed  burst  error-correcting  capability  of  RS  codes  is  first  explored  for  sin¬ 
gle  and  multiple  burst  errors.  This  burst  error-correcting  capability  is  then  used  to  define 
the  burst  error-correcting  efficiency  of  selected  RS  codes.  The  idea  of  guaranteed  burst 
error  correction  and  burst  error  correction  efficiency  is  then  expanded  to  include  block 
and  convolutional  interleaving  techniques.  Many  useful  tables  and  graphs  of  guaranteed 
burst  error-correcting  capabilities  and  efficiencies  have  been  generated  for  the  code  de¬ 
signer. 
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2  Introduction 


Why  the  concern  for  burst  error  correction?  Other  than  deep  space  communications, 
an  overwhelming  majority  of  practical  communication  systems  cannot  be  modeled  as 
simple  Additive  White  Gaussian  Noise  (AWGN)  channels.  In  recording  systems, 
whether  they  are  magnetic  or  optical,  defects  in  the  recording  material  and  recording 
mechanisms  cause  systematic  errors.  Wireless  communication  systems  with  fading  ef¬ 
fects  are  also  faced  with  the  occurrence  of  non-random  errors.  Jamming  environments, 
especially  intentional  jamming  as  in  military  applications,  are  confronted  with  the  effects 
of  non-random  errors.  Communication  channels  that  have  non-random  characteristics  are 
considered  to  contain  memory  and  tend  to  produce  what  is  known  as  burst  errors.  In  fact, 
most  operational  communication  system  applications  exhibit  a  combination  of  random 
and  non-random  characteristics.  These  channels  are  termed  diffuse  or  compound  chan¬ 
nels  [Swe]  and  cannot  be  modeled  as  AWGN  channels.  Therefore,  in  answering  the 
question  why  burst  error  correction,  most  real  world  communication  systems  need  codes 
that  can  combat  the  common  occurrence  of  burst  errors. 
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2.1  Burst  Errors 


The  purpose  of  this  section  is  to  provide  the  reader  with  a  better  understanding  of  the 
errors  associated  with  digital  communication  systems.  First  a  simple  example  and  defi¬ 
nition  of  errors  that  occur  in  binary  codewords  is  given.  A  discussion  of  the  two  different 
types  of  possible  errors  follows.  Then  the  definition  of  errors  is  expanded  to  include  the 
word  oriented  RS  codes.  The  section  is  concluded  with  brief  examples  of  the  classic  er¬ 
ror  causing  bursty  channel. 

A  received  encoded  block  of  data  can  be  interpreted  as  the  linear  combination  of  a 
valid  codeword  and  an  error  pattern.  This  relationship  is  shown  as  R  =  C  +  E  where  R  is 
equal  to  the  received  block  of  data,  C  is  a  valid  codeword  and  E  is  an  error  pattern  and  the 
summation  implies  binary  vector  addition.  In  the  binary  case  a  ‘  1  ’  in  the  error  pattern 
represents  the  occurrence  of  an  error  in  that  position.  As  an  example,  consider  a  commu¬ 
nication  system  where  the  sixteen  bit  block  of  received  data  is  R  =  0010001 110010000. 
Then  suppose  the  original  transmitted  data  was  the  valid  codeword  C  = 
0001010010001000.  Therefore,  because  of  errors  in  the  channel,  an  error  pattern  E  = 

001 101 1 10001 1000  is  produced.  It  can  be  observed  that  the  modulo  two  addition  of  R 
and  E  will  correct  the  errors  in  R  and  return  the  intended  block  of  data  C. 

Errors  occur  in  communication  systems  in  two  different  forms,  random  and  burst  er¬ 
rors.  An  example  of  a  random  error  is  given  in  Figure  2.0  where  n  bits  of  a  received  data 
stream  are  given  and  x  represents  a  bit  error.  There  are  five  random  errors  shown  in  Fig¬ 
ure  2.0  and  the  error  pattern  would  be  E  =  0001010100011000  with  weight  e  equal  to 
five.  Note  that  for  no  known  deterministic  reason,  two  of  the  random  errors  occurred  ad- 
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jacent  to  each  other.  Additive  white  Gaussian  noise  channels  are  random  error  channels. 


Figure  2.0  Random  Errors 

Burst  errors  are  conventionally  defined  in  terms  of  their  length,  and  can  be  viewed  as  a 
cluster  of  individual  errors.  The  length  of  a  burst  error  is  equal  to  the  number  of  bits  that 
span  the  error  cluster.  Therefore,  a  burst  error  of  length  b  is  defined  as  the  number  of  bits 
in  the  error  cluster  starting  from  the  first  bit  in  error  (first  ‘  1’  in  the  error  pattern)  and 
ending  with  the  last  bit  in  error  (last  ‘  1’  in  the  error  pattern)  [Lin],  For  example  the  burst 
error  shown  in  Figure  2.1  is  of  length  b  =  1 1  bits  and  has  an  error  pattern  of  E  = 

0011011 10001 1000  with  weight  e  equal  to  seven. 


Figure  2.1  Conventional  Burst  Error 

More  than  one  burst  error  may  be  imposed  onto  a  codeword.  When  this  happens,  the 
codeword  is  said  to  contain  multiple  burst  errors.  For  example  in  Figure  2.1  above,  the 
codeword  may  be  considered  to  contain  three  burst  errors  of  which  two  are  of  weight  two 
and  one  burst  error  of  weight  three.  In  keeping  with  the  notation  used  in  [Ben],  a  code¬ 
word  may  contain  s  burst  errors  (s  >  1)  of  length  b  bits  (b  >  1)  and  each  burst  of  weight  e 
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(e>  1). 


To  resolve  any  ambiguities  about  defining  burst  errors,  the  idea  of  a  guard  space  is 
now  explained.  Just  as  a  channel  with  memory  produces  error  clusters,  the  channel  also 
produces  spaces  where  no  errors  occur.  These  clusters  of  error  free  bits  are  indicated  as 
sequences  of  all  zeros  in  the  error  pattern  E  and  are  called  guard  spaces.  Therefore,  a 
burst  error  b  is  defined  in  relation  to  the  surrounding  guard  space  g  [Gall].  For  example, 
in  Figure  2.2  the  knowledge  of  the  guard  spaces  dictates  the  definition  of  the  possible 
burst  errors.  The  bits  of  information  indicated  by  x’s  in  Figure  2.2  can  represent  burst 
errors  of  length  less  than  or  equal  to  b. 


XXX  x  |  000  . 0 

XXX.  x  |  000.  .  0 

XXX  .  .X 

000  0 

Figure  2.2  Burst  Errors  with  Guard  Spaces 

The  definition  of  burst  errors  will  now  be  expanded  for  the  case  of  RS  codes  that  are 
examined  in  this  report.  Reed-Solomon  codes  are  word-oriented  codes,  in  which  every 
codeword  is  composed  of  m  bit  symbols.  We  assume  that  the  random  error-correcting 
capability  of  a  RS  code  is  denoted  by  t.  If  any  one  of  the  m  bits  of  a  given  symbol  is  in 
error,  then  the  entire  symbol  is  considered  to  be  erroneous.  Therefore,  it  is  customary  to 
measure  random  errors  in  a  RS  code  according  to  the  number  of  symbols  in  error  and  not 
the  number  of  bits.  Figure  2.3,  a  representation  of  a  n  symbol  codeword,  indicates  two  (t 
=  2)  random  errors. 
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Figure  2.3  Word  Oriented  Random  Error 


Let’s  assume  a  burst  error  of  length  b  bits  corrupts  an  RS  code  and  t  consecutive  sym¬ 
bol  errors  occur.  To  define  the  relation  between  b  and  t  for  burst  errors  in  RS  codes,  as¬ 
sume  all  burst  errors  are  of  the  worst  case.  Here  the  worst  case  burst  error  is  defined  as  a 
burst  error  in  which  all  b  bits  of  the  error  cluster  are  in  error.  The  worst  case  burst  error 
would  have  an  error  pattern  with  a  burst  pattern  of  all  ‘  l’s  (e.g.  E  =  00001 111111 10000 
has  burst  error  length  of  b  =  8  bits  as  shown  in  Figure  2.4).  If  the  burst  error  occurs  ex¬ 
actly  over  t  symbols  as  shown  in  Figure  2.4,  then  the  length  of  the  burst  is  equal  to  the 
number  of  bits  over  t  symbols  (b  =  t-m  bits).  However,  a  burst  error  of  the  same  length 
can  corrupt  up  to  t  +  1  symbols  as  shown  in  Figure  2.5.  The  burst  error-correcting  capa¬ 
bility  of  RS  codes  will  be  further  examined  in  Section  4  of  this  document. 
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Figure  2.4  Word  Oriented  Burst  Error  Located  Over  Exactly  ( t )  Symbols 
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Figure  2.5  Word  Oriented  Burst  Error  Over  (t  +  1)  Symbols 


2.2  The  Bursty  Channel 

There  are  several  systems  of  interest  that  operate  over  the  classic  bursty  channel  which 
use  burst  error-correcting  RS  codes.  Burst  errors  are  found  in  both  magnetic  and  optical 
recording  systems.  In  digital  mobile  communication  systems,  such  as  cellular  and  satel¬ 
lite  communication  systems,  burst  errors  are  a  major  concern.  And  in  military  communi¬ 
cations,  intentional  jamming  may  introduce  burst  errors  into  the  channel.  Burst  errors  can 
have  serious  adverse  effects  on  signal  reliability,  and  hence  the  need  for  burst  error-cor¬ 
recting  codes. 

With  increases  in  recording  speeds  and  data  storage  densities  for  magnetic  recording 
systems,  burst  errors  have  become  common  place.  Dust  and  defects  in  the  recording  ma¬ 
terial  and  the  read/write  heads  of  a  magnetic  disk  recording  system  are  usually  to  blame 
for  the  occurrence  of  burst  errors.  These  defects  have  been  known  to  produce  burst  errors 
of  hundreds  of  bits  in  length  [Lin,  chap  16].  One  example  of  the  use  of  RS  codes  in  mag¬ 
netic  recording  systems  is  the  IBM  3370  disk  storage  system  [Mee]. 

Optical  recording  systems  have  used  RS  codes  extensively  for  the  correction  of  burst 
errors.  One  digital  recording  system  that  owes  its  existence  to  RS  codes  is  the  compact 
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disk  (CD)  recording  system.  With  the  reduction  in  size  of  the  recording  media  (a  com¬ 
pact  disk  has  an  outer  diameter  of  120mm  [Car]),  the  density  of  the  data  has  become 
enormous.  The  track  width  of  a  CD  is  roughly  0.6  pm  and  one  bit  of  audio  data  covers 
1 .0  pm2  of  disk  space  [Car].  In  comparison  to  the  LP  record  track,  the  CD  has  a  track 
ratio  of  approximately  60:1  [Bae].  That  is,  sixty  tracks  can  be  placed  on  a  CD  in  the 
same  amount  of  space  once  used  by  one  LP  track.  With  these  data  densities,  it  is  clear 
that  the  slightest  defect  or  scratch  may  cause  extremely  long  burst.  If  it  were  not  for  the 
error  control  scheme  known  as  Cross  Interleaved  Reed-Solomon  Code  (CIRC),  the  Com¬ 
pact  Disk  recording  system  could  not  deliver  nearly  the  quality  of  music  it  now  does. 

The  CIRC  error  correction  scheme  has  the  capability  to  correct  burst  errors  of  approxi¬ 
mately  4000  bits  in  length,  this  corresponds  to  roughly  2.5  mm  of  track  length  [Imm2]. 
For  a  complete  description  of  the  CIRC  error  control  method  as  applied  to  the  compact 
disk  refer  to  [Vri]  and  [Imml]. 

Since  the  onset  of  the  digital  communication  era,  mobile  communication  has  become  a 
fast  developing  area.  Popular  topics  in  the  mobile  communication  arena  include  cellular 
telephone  systems,  personal  communication  systems  (PCS),  and  satellite  communication 
systems.  All  of  these  communication  systems  are  faced  with  channel  problems  such  as 
multipath  effects,  channel  fading  and  Doppler  shifts.  In  light  of  such  adverse  channel 
impairments,  burst  errors  are  very  common.  One  example  of  the  use  of  RS  codes  in  mo¬ 
bile  communications  is  in  the  cellular  digital  packet  data  (CDPD)  system  [Kami] 

[Kam4],  For  a  more  detailed  review  of  the  signal  characteristics  of  a  mobile  communi¬ 
cations  channel,  see  [Vau]. 
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Fading  in  a  mobile  communication  channel  cam  be  statistically  modeled  as  a  Rayleigh 
distributed  process  [Wic3].  Fades  are  categorized  into  what  are  known  as  slow  and  fast 
fades.  The  size  of  the  burst  error  produced  is  directly  related  to  the  length  of  the  fade  im¬ 
posed  on  the  system.  Slow  fades  are  longer  in  duration  relative  to  the  symbol  rate  of  the 
communication  signal  being  affected.  Therefore,  a  slow  fade  will  produce  multiple  sym¬ 
bol  errors.  It  has  been  shown  that  mobile  communication  systems  with  the  aid  of  error 
control  coding  and  interleaving  techniques  have  improved  performance  [Yue]. 

In  today’s  military  combat  environment,  the  ability  to  communicate  determines  win¬ 
ners  or  losers.  For  this  reason  electronic  warfare  (EW)  and  information  warfare  (IW)  are 
targets  of  many  discussions  and  much  research.  Of  particular  interest  is  the  survival  of  a 
communication  system  in  the  presence  of  intentional  jamming.  There  are  many  different 
types  of  jamming  techniques,  however  pulsed  jamming  has  shown  to  be  the  most  optimal 
for  creating  interference  [Tor].  Pulsed  jamming  concentrates  more  of  the  total  signal 
power  into  each  pulse  unlike  a  continuous  wave  jammer  whose  power  is  more  spread  out. 
If  the  pulse  width  of  the  jamming  signal  is  of  the  same  order  as  the  symbol  transmission 
rate,  then  approximately  one  symbol  will  be  affected  per  jamming  pulse.  If  the  pulse 
width  is  greater  than  the  symbol  transmission  rate,  then  multiple  symbols  can  be  affected 
per  jamming  pulse  causing  burst  errors.  As  illustrated  in  Figure  2.6  below,  jamming  sig¬ 
nal  J,(t)  has  pulse  width  approximately  equal  to  the  symbol  rate  of  the  communication 
signal  S(t).  Jamming  signal  J2(t)  has  a  pulse  width  approximately  equal  to  three  times 
the  symbol  rate  of  S(t)  thus  causing  burst  errors. 
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S(t) 


Figure  2.6  Pulsed  Jamming 

Techniques  such  as  spread  spectrum  communication  have  been  employed  with  good 
success  in  defending  against  intentional  jamming.  However,  it  can  be  shown  that  using 
spread  spectrum  with  the  addition  of  error  control  coding  and  interleaving  techniques,  the 
affect  of  pulsed  noise  jamming  can  almost  be  eliminated  [Zie,  chap  1 1&12]  [Pur].  One 
example  of  a  military  communication  system  that  uses  RS  codes  is  the  Joint  Tactical  In¬ 
formation  Distribution  System  (JTIDS)  being  supported  by  all  branches  of  the  military 
and  NATO  [Zie,  chap  13]. 

With  the  ever  increasing  complexity  of  communication  systems,  the  powerful  burst 
error-correcting  RS  codes  will  no  doubt  see  more  applications  in  digital  recording  sys¬ 
tems,  mobile  communication  systems,  and  in  communication  environments  with  interfer¬ 
ence. 
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3  Finite  Field  Mathematics 


To  understand  Reed-Solomon  codes  an  introduction  to  Galois  (pronounced  ‘gal-wa’) 
field  algebra  is  presented.  To  learn  and  understand  Galois  fields  it  is  important  to  have  a 
foundation  in  the  laws  of  abstract  algebra  that  govern  Galois  field  theory.  However,  it  is 
by  no  means  the  intention  of  this  section  to  present  an  in-depth  coverage  of  abstract  alge¬ 
bra.  The  purpose  of  this  section  is  to  present  an  overview  of  basic  ideas  and  laws  needed 
in  understanding  Galois  field  algebra  as  applicable  to  RS  coding.  The  reader  is  encour¬ 
aged  to  review  any  of  the  following  references  for  a  more  detailed  analysis  of  abstract 
algebra  and  Galois  field  theory  [Cla],  [Lin],  [McEl],  [Petl],  and  [Wicl]. 

3.1  Groups 

A  group  (G)  is  a  set  in  which  a  binary  operation  ( *  )  has  been  defined.  Binary  in  the 
sense  that  the  operation  is  performed  on  two  elements.  The  resultant  of  the  binary  opera¬ 
tion  is  also  another  element  in  the  set.  This  unique  property  is  called  closure.  Three 
other  properties  are  required  in  order  for  a  set  to  form  a  group.  The  law  of  associativity 
must  be  upheld.  There  must  exist  an  identity  element  in  the  group.  Every  element  in  the 
group  must  have  a  unique  inverse.  Optionally,  if  a  group  satisfies  as  a  fifth  property  the 
law  of  commutativity,  the  group  is  said  to  be  commutative  and  called  an  Abelian  group. 
In  summary,  the  follow  properties  apply: 

a)  Closure:  V  a,b  e  G;  a  •  b  =  c;  where  c  e  G 

b)  Associativity:  (a  •  b)  •  c  =  a  •  (b  •  c) 

c)  Identity:  V  a  e  G  there  exist  e  e  G  such  that  a  •  e  =  e  •  a  =  a 
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d)  Inverse:  V  a  e  G  there  exist  a'1  e  G  such  that  a  •  a'1  -  a'1  •  a  -  e 

e)  Commutativity:  V  a,b  e  G;  a  •  b  =  b  •  a  (optional) 

3.2  Rings 

A  ring  ( R )  is  a  set  of  elements  with  two  binary  operations  ( +  and  • )  defined.  There 
are  three  properties  that  must  hold  for  a  set  to  form  a  ring.  The  ring  must  form  an  Abe¬ 
lian  group  under  ( + )  with  the  additive  identity  element  indicated  by  ( 0 ).  The  law  of 
associativitity  must  be  upheld  for  the  binary  operation  ( • ).  Also,  the  binary  operation  of 
( • )  must  be  distributive  over  the  binary  operation  of  ( + ).  Additionally,  if  the  operation 
(  •  )  is  commutative,  then  the  ring  is  called  a  commutative  ring.  If  the  operation  ( •  )  has 
an  identity  element  labeled  as  (  1 ),  then  the  ring  is  said  to  be  a  ring  with  identity.  If  the 
ring  meets  both  of  the  additional  properties,  then  it  is  called  a  commutative  ring  with 
identity.  The  ring  is  an  important  building  block  for  error  control  coding.  It  is  the  ring 
that  allows  for  the  development  of  the  polynomial.  Summarizing,  the  following  proper¬ 
ties  apply: 

a)  Commutative  group  under  addition  (  +  )  with  identity  element  (  0 ) 

b)  Associativity  under  multiplication  (  •  ):  V  a,b,c  e  R;  (a  •  b)  *  c  =  a  •  (b  •  c) 

c)  Multiplication  (  •  )  distributes  over  addition  (  +  ):  V  a,b,c  e  R; 
a  •  (b  +  c)  =  (a  •  b)  +  (a  •  c) 

d)  Commutativity  under  multiplication  ( •  ):  V  a,b  e  R;  a  •  b  =  b  •  a 

e)  Multiplication  ( *  )  has  an  identity  element  labeled  as  ( 1  ) 

3.3  Fields 
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Now  using  the  concept  of  the  group,  the  algebraic  structure  of  a  field  can  be  intro¬ 
duced.  A  field  ( F )  is  defined  such  that  it  can  be  thought  of  as  a  combination  of  a  com¬ 
mutative  group  under  addition  ( +  )  and  a  commutative  group  under  multiplication  (  • ). 
Basically,  this  makes  the  field  a  set  of  elements  in  which  the  operations  of  addition  (  +  ), 
multiplication  ( •  ),  subtraction  ( - ),  and  division  ( -r  )  are  defined.  The  follow  properties 
must  be  satisfied  to  define  a  field: 

a)  F  is  a  commutative  group  under  ( + )  with  additive  identity  element  (  0 ). 

b)  The  set  of  nonzero  elements  in  F  is  a  commutative  group  under  ( • ) 
with  multiplicative  identity  element  (  1  ). 

c)  Multiplication  ( • )  distributes  over  addition  (  +  ),  that  is  V  a,b,c  €  F; 
a  •  (b  +  c)  =  (a  •  b)  +  (a  •  c) 

In  block  coding,  fields  of  finite  order  are  key  algebraic  structures  of  interest.  A  finite 
field  is  referred  to  as  a  Galois  field  that  is  named  after  Evariste  Galois,  the  discoverer  of 
finite  fields.  A  Galois  field  is  denoted  by  GF(q),  where  the  order  of  the  field  is  q. 

It  should  be  obvious  from  the  definition  that  a  field  needs  to  contain  at  least  two  ele¬ 
ments,  the  additive  identity  element  ( 0  )  and  the  multiplicative  identity  element  ( 1  ). 

The  field  that  is  formed  from  these  two  elements  is  known  as  the  binary  field  GF(2).  The 
binary  field  and  its  extensions  are  very  important  and  have  many  applications  in  engi¬ 
neering.  Before  presenting  the  binary  extension  field,  primitive  elements  and  primitive 
polynomials  must  be  introduced.  But  first  the  concept  of  a  vector  space  will  be  intro¬ 
duced. 
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3.4  Vector  Spaces 


Let  V  represent  a  set  of  elements  called  vectors  and  F  represent  a  field  of  elements 
called  scalars.  Two  additional  operations  are  introduced  that  can  be  applied  to  the  ele¬ 
ments  of  V  and  F.  Let  ( +  )  be  a  binary  additive  operation  that  maps  linear  combinations 
of  vectors  v„  v2,  v3, contained  in  V  into  a  single  vector  Vj  also  contained  in  V.  This 
binary  additive  operation  is  referred  to  as  vector  addition.  And  let  ( • )  represent  a  binary 
multiplicative  operation  performed  on  a  scalar  a,  contained  in  F  and  a  vector  Vj  contained 
in  V,  such  that  ( * )  maps  a,  and  Vj  into  a  vector  u;  =  a;  •  vs  also  contain  in  V.  A  vector 
space  is  formed  if  the  following  rules  are  meet: 

a)  V  forms  a  commutative  group  under  ( + ). 

b)  V  a  e  F  and  VveV;  a  •  v  =  u  e  V. 

c)  Associativity  under  multiplication  ( •  ):  V  a,b  e  F  and  VveV; 

(a  •  b)  *  v  =  a  *  (b  •  v) 

d)  Multiplication  (  * )  distributes  over  addition  (  +  ):  V  a,b  e  F  and  V  u,v  e  F; 
a  •  (u  +  v)  =  (a  •  u)  +  (a  •  v),  and  (a  +  b)  •  u  =  (a  •  u)  +  (b  •  u). 

e)  V  v  e  V,  1  •  v  =  v  •  1  =  v;  the  scalar  multiplication  identity  element  is  1. 

In  a  vector  space,  F  is  commonly  referred  to  as  the  “ground  field”  or  “base  field”  of 
the  vector  space.  The  properties  of  the  vector  space  allow  for  the  existence  of  what  is 
known  as  a  vector  n-tuple.  A  n-tuple  has  the  form  v  =  (v0,  v„  v2,  ...,vn.,)  where  the  ele¬ 
ments  (vj)  are  from  the  ground  field  F.  It  is  the  n-tuple  representation  of  a  vector  space 
that  is  so  important  in  implementing  error  control  codes  due  to  the  association  of  binary 
words  with  vector  n-tuples.  For  example,  if  the  ground  field  is  taken  to  be  binary 
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(  GF(2) ),  then  00101  is  a  five  bit  binary  word  can  be  represented  by  v  =  (00101)  a  5- 
tuple  taken  from  the  vector  space  V5 .  Of  course,  here  the  vector  space  is  finite  and  com¬ 
posed  of  32  (25)  vector  5-tuples. 

Some  other  important  ideas  of  vector  spaces  associated  with  error  control  coding  are 
the  dimension  of  the  vector  space,  the  basis  of  a  vector  space,  and  the  spanning  set  of  a 
vector  space.  A  spanning  set  S  =  (v0,  v„  v2,  ...,vn.,),  is  a  set  of  vectors  from  which  the  lin¬ 
ear  combination  of  the  vectors  (ao*v0  +  a,’v,  +  a2*v2  +...+  an.,*vn.l ;  where  a,  =  0  or  1) 
forms  all  other  possible  vectors  of  the  vector  space.  Going  another  step  further,  a  basis  is 
a  set  of  vectors  of  minimum  order  that  spans  the  vector  space.  Therefore,  a  basis  is  a 
spanning  set  that  has  no  additional  redundant  vectors.  A  basis  set  is  always  a  spanning 
set  but  a  spanning  set  may  not  be  a  basis.  Finally,  the  dimension  of  a  vector  space  is 
equal  to  the  number  (k)  of  vectors  in  the  basis  set.  This  of  course  leads  to  the  fact  that 
each  vector  of  a  vector  space  is  unique. 

3.5  Primitive  Element 

Elements  of  a  finite  field  have  individual  cardinalities.  Defining  a2  =  a  •  a,  a3  =  a  •  a  • 
a,  etc...,  the  order  of  an  element  (  n  =  ord(a) )  is  the  value  that  satisfies  a"  =  1 ,  the  identity 
element.  Of  particular  importance,  are  the  elements  in  field  GF(q)  with  order  q-1,  which 
are  called  primitive  elements.  It  can  be  proven  that,  since  a  primitive  element  has  order 
equal  to  one  less  than  that  of  the  size  of  the  field,  consecutive  powers  of  a  primitive  ele¬ 
ment  can  represent  all  nonzero  elements  of  a  field.  It  also  has  been  proven  that  there  exist 
at  least  one  element  in  GF(q)  with  order  q-1,  and  hence  there  is  always  at  least  one  primi- 
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tive  element  associated  with  a  finite  field.  The  importance  of  primitive  elements  will  be 
further  illustrated  in  the  section  on  the  development  of  binary  extension  fields. 

3.6  Primitive  Polynomials 

It  has  been  shown  thus  far  that  a  field  GF(q)  can  be  represented  by  0  and  q-1  consecu¬ 
tive  powers  of  a  primitive  element  (  a  ).  Remembering  the  additive  and  multiplicative 
nature  of  the  field,  a  field  can  be  composed  of  polynomials  of  mlh  degree  whose  coeffi¬ 
cients  are  taken  from  a  defined  finite  field.  The  notation  GF(q)[x]  represents  such  a 
polynomial  f(x)  with  coefficients  taken  from  GF(q).  As  defined  earlier,  this  collection  of 
polynomials  is  called  a  commutative  ring  with  identity. 

A  polynomial  f(x)  which  cannot  be  factored  into  the  product  of  lower  degree  polyno¬ 
mials  in  GF(q)[x],  is  an  irreducible  polynomial  in  GF(q).  That  is  to  say,  the  polynomial 
f(x)  is  in  the  simplest  form  over  GF(q)[x].  For  example,  f(x)  =  x2  +  x  +  1  is  irreducible  in 
GF(2)[x],  however  f(x)  can  be  factored  in  GF(4)[x].  With  this  illustration,  it  is  clear  that 
when  defining  irreducibility  the  field  must  be  specified. 

An  irreducible  polynomial  p(x)  of  arbitrary  degree  m,  defined  in  GF(p)[x]  where  p  is  a 
prime,  is  called  a  primitive  polynomial  if  the  smallest  value  of  n  for  which  p(x)  divides 
xn  -  1  is  n  =  pm  - 1 .  For  example,  p(x)  =  x3  +  x  +  1  is  a  primitive  polynomial  in  GF(2)[x] 
and  the  smallest  degree  polynomial  of  the  form  x"  - 1  that  p(x)  divides  is  x7  - 1 .  Note 
that  all  primitive  polynomials  in  GF(q)[x]  are  irreducible,  but  not  all  irreducible  polyno¬ 
mials  are  primitive.  An  extremely  important  property  of  primitive  polynomials,  which 
will  be  exploited  in  the  development  of  RS  codes,  is  that  the  roots  (a1)  of  a  primitive 
polynomial  p(x)  in  GF(p)[x]  have  order  pm  -  1  .  Recall,  an  element  of  order  pm  - 1 
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forms  all  nonzero  elements  of  a  field  GF(q)  of  order  q  =  pm  - 1 .  Summarizing,  an  mth  de¬ 
gree  polynomial  in  GF(p)[x],  is  a  function  of  x  with  coefficients  in  GF(p)  and  has  roots 
(a1)  over  the  field  GF(pm  ). 

3.7  Binary  Extension  Fields 

In  the  study  of  error  control  coding  we  are  interested  in  polynomials  in  GF(p)[x]  that 
have  coefficients  in  GF(p)  and  roots  (a1)  in  GF(pm),  where  p  is  equal  to  the  binary  ground 
field.  Simply  stated,  p  =  2.  By  defining  a  field  in  this  fashion,  the  elements  of  the  field 
become  binary  n-tuples  of  length  m  over  GF(pm).  Essentially,  binary  words  of  length  m 
are  defined.  Using  a  primitive  polynomial  p(x)  of  degree  m,  and  the  fact  a  is  a  primitive 
root  of  p(x),  the  construction  of  a  binary  extension  field  of  order  pm  will  now  be  illus¬ 
trated. 

a)  Step  1 :  Choose  a  primitive  polynomial  of  degree  m,  where  m  is  equal  to  the 
size  of  the  binary  word  of  interest.  For  instance,  if  working  with  3-bit  binary 
symbols,  choose  a  primitive  polynomial  of  degree  three.  Choose  p(x)  =  x3  +  x 
+1. 

b)  Step  2:  Let  a  be  a  primitive  root  of  p(x).  Therefore,  p(a)  =  a3  +  a  +1  =  0  and 

a3  =  a  +  1 . 

c)  Step  3:  Recall  that  p(x)  must  divide  f(x)  =  x7  -  1 ;  where  n  =  23  -  1 .  If  p(x)  di¬ 
vides  f(x)  and  has  a  as  a  root,  then  a  must  also  be  a  root  of  f(x).  Therefore, 

f(a)  =  a7  - 1  =  0  and  hence  a7  =  1 . 
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d)  Step  4:  Include  the  zero  element.  Since  the  consecutive  powers  of  the  primi¬ 
tive  element  a  only  form  the  nonzero  elements  of  the  field,  the  zero  element 
must  be  included. 

e)  Step  5:  Begin  construction  of  the  binary  extension  field  using  the  properties 
of  the  above  steps  to  keep  all  elements  in  their  reduced  format.  The  field  ele¬ 
ments  are 

0,  a°=l,  a‘=a,  a2,  a3  =  a+l,  a4  =  a3a  =  (a+1) a  =  a2  +  a, 
a5  =  a4a  =  (a2  +  a)a  =  (a3  +  a2)  =  (a  +  1  +  a2)  =  a2  +  a  +  1,  etc... 


Table  3.0  Construction  of  Binary  Extension  Field  GF(8) 


EXPONENTIAL  POLYNOMIAL  BINARY  N-TUPLE 

REPRESENTATION  REPRESENTATION  REPRESENTATION 


0 

0 

000 

1 

1 

100 

a 

a 

010 

a2 

a2 

001 

a3 

a  +  1 

no 

a4 

a2  +  a 

Oil 

a5 

a2  +  a  +  1 

in 

a6 

a2  +  1 

101 
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4  Reed-Solomon  Codes 


With  the  evolution  of  today’s  complex  communication  systems,  it  is  apparent  that  a 
strong  error-correcting  scheme  is  needed  to  combat  burst  errors.  Reed-Solomon  codes, 
once  thought  of  as  only  a  great  non-practical  mathematical  theory,  have  proven  to  be  the 
answer  to  the  burst  error-correcting  problem.  Reed-Solomon  codes  received  their  glory 
from  such  applications  as  the  deep  space  exploration  missions  of  Galileo  and  Voyager 
[McE2],  the  compact  disc  optical  recording  system  [Imm2],  and  many  applications  in 
magnetic  recording  systems  [Her].  With  the  advances  being  made  in  communication 
systems,  RS  codes  are  sure  to  see  many  future  applications.  The  remainder  of  this  section 
will  attempt  to  explain  what  RS  codes  are  and  how  they  are  encoded  and  decoded. 

Reed-Solomon  codes  are  a  subclass  of  linear  block  codes.  Like  any  other  block  codes, 
the  RS  codes  append  parity  information  to  a  block  of  data  so  that  at  the  receiver  error  de¬ 
tection  and  correction  is  possible.  Reed-Solomon  codes  are  typically  specified  by  their 
length  (n)  and  dimension  (k).  Therefore,  we  specify  a  particular  RS  code  as  a  (n,k)  RS 
code.  Recall  from  Section  3  that  k  is  the  number  of  elements  in  the  basis  of  the  vector 
space  and  is  an  indication  of  the  number  of  data  symbols  in  the  original  information 
block.  If  n  (where  n  >  k)  is  equal  to  the  number  of  symbols  in  the  encoded  block,  then 
there  is  n-k  parity  symbols  appended  to  the  k  original  data  symbols.  Also,  a  RS  code  is 
said  to  have  a  code  rate  defined  as  Rc  =  k  /  n.  Refer  to  Figure  4.0  below  describing  the 
structure  of  a  (7,3)  RS  codeword. 
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«- 


n  symbols 


m  bits — > 


i - n  -  k  parity  symbols - ►  «■ - k  data  symbols - ► 

Figure  4.0  Structure  of  a  (7,3)  Reed-Solomon  Code 

4.1  Error-correcting  Capability 

Reed-Solomon  codes  are  word  oriented  linear  block  codes  capable  of  correcting  t  ran¬ 
dom  symbol  errors.  Reed-Solomon  codes  have  the  ability  to  correct  burst  errors  of  length 
equal  to  t  consecutive  symbols.  Since  RS  codes  are  word  oriented,  that  is  they  are  com¬ 
posed  of  m-bit  symbols,  they  can  correct  bursts  of  maximum  length  equal  to  t  •  m  bits.  It 
can  be  shown  that  the  guaranteed  length  for  burst  error  correction  of  a  RS  code  is 
b  =  m(t  -  1)  +  1  bits.  This  can  be  expanded  to  guaranteed  double  burst  error  correction 
(b  =  m(Lt/2j  - 1)  +  1)  and  guaranteed  triple  burst  error  correction  (b  =  m(Lt/3j  - 1)  +  1) 
[Kam2].  The  definitions  given  here  for  the  guaranteed  (not  maximum)  burst  error  cor¬ 
rection  will  be  utilized  in  the  following  sections  to  define  the  efficiencies  of  RS  codes. 

4.2  Maximum  Distance  Codes 

Reed-Solomon  codes  have  the  important  property  of  being  Maximum  Distance  Sepa¬ 
rable  (MDS)  codes.  It  has  been  shown  that  RS  codes  satisfy  the  Singleton  bound  with 
equality  and  therefore  are  considered  Maximum  Distance  codes  [Wicl].  Maximum  dis- 
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tance  codes  have  the  highest  obtainable  error  correction  capability  of  any  block  code  of 
the  same  redundancy  [Sin]. 

Hamming  weight  has  been  defined  as  the  number  of  non-zero  elements  in  a  binary 
symbol.  For  example,  the  symbol  (1 1010001)  has  a  Hamming  weight  of  four  (denoted  W 
=  4).  Hamming  distance  is  defined  as  the  number  of  positions  between  two  binary  sym¬ 
bols  that  are  in  difference.  For  example,  the  Hamming  distance  between  the  two  binary 
symbols  (11010001)  and  (10001010)  is  equal  to  five  (denoted  D  =  5).  The  Hamming 
distance  is  also  shown  to  be  the  Hamming  weight  of  the  modulo-two  addition  of  any 
combination  of  symbols.  For  example,  (11010001)  +  (10001010)  =  (0101 1011)  and  the 
Hamming  weight  of  the  resultant  (01011011)  is  five  as  before.  The  minimum  distance 
(<!„,;„)  of  a  code  is  defined  as  the  minimum  Hamming  distance  between  all  possible  valid 
pairs  of  codewords  in  a  given  vector  space.  Recall  that  a  vector  space  defined  over  a  bi¬ 
nary  extension  field  is  composed  2"  possible  code- vectors,  where  n  equals  the  code 
length.  For  the  same  vector  space,  there  are  2k  valid  codewords,  where  k  is  the  dimension 
of  the  vector  space.  Therefore,  there  are  at  least  dmin  bits  separating  valid  codewords. 
Maximum  distance  codes,  are  codes  that  have  the  largest  possible  minimum  distance  of 
separation  between  valid  codewords.  In  fact  it  has  been  shown  by  [Sin]  that  if  a  code  has 
d^  =  n  -  k  +  1 ,  it  is  a  MDS  code.  This  can  be  explained  by  visualization  of  any  two 
codewords  n  symbols  in  length  with  k  data  symbols  and  n-k  parity  symbols.  To  be  valid 
codewords,  the  k  data  symbols  must  at  minimum  differ  in  one  position.  Also,  the  most 
desired  case  would  be  when  all  n-k  parity  symbols  differ  between  two  valid  codewords. 
Combining  the  above  two  statements,  the  largest  minimum  separation  between  valid 
codewords  is  equal  to  the  maximum  possible  distance  between  all  parity  symbols 
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(D  =  n-k)  and  the  single  distance  between  data  symbols  ( D  =  1),  yields  =  n  -  k  +  1 . 


4.3  Encoding 

Since  the  introduction  of  RS  codes  in  1960  by  Irving  Reed  and  Gus  Solomon  [Ree], 
there  have  been  many  different  approaches  for  encoding  of  RS  codes.  These  differences 
are  due  to  the  fact  that  there  existed  no  efficient  methods  of  decoding  RS  codes,  therefore 
new  encoding  methods  where  created  that  gave  way  to  better  decoding  algorithms.  The 
original  method  defining  RS  codes  as  polynomials  over  certain  finite  fields  can  be  re¬ 
viewed  in  Reed  and  Solomon’s  original  1960  paper  [Ree].  There  is  also  a  summary  in 
[Wic2].  The  method  of  encoding  RS  codes  that  will  follow  is  the  extended  non-binary 
BCH  code  approach.  There  is  a  reason  for  choosing  the  extended  BCH  approach.  In 
most  applications  this  approach  leads  to  more  efficient  decoding  algorithms. 

In  order  to  construct  an  RS  code  as  an  extension  of  a  BCH  code,  one  must  use  the 
generator  polynomial  approach.  RS  codes  are  represented  by  generator  polynomials  g(x). 
These  polynomials  are  constructed  over  GF(2m)  and  therefore  have  as  roots  the  primitive 
elements  a' .  In  order  to  be  MDS,  g(x)  must  meet  the  BCH  bound.  Therefore,  the  a1 
primitive  elements  must  be  2t  =  n-k  consecutive  primitive  roots  of  g(x).  Such  a  polyno¬ 
mial  must  be  of  the  form  g(x)  =  (x  -  ab )  (x  -  ab+1 )  (x  -  a1*2 )...  (x  -  ab+2t'] ).  For  narrow- 
sense  RS  codes,  g(x)  can  be  represented  by  the  following  equation  (in  which  b  equals 
one): 

2t 

g(x)  =  n  (x  -  a' ) 
i=l 
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When  all  the  factors  of  g(x)  are  multiplied  out  to  form  the  generator  polynomial,  it  will  be 
a  polynomial  of  degree  2t  =  n  -  k  defined  over  the  field  GF(q)  where  q  =  2m  and  will  thus 
have  qk  valid  codewords. 

By  defining  RS  codes  as  extended  BCH  codes  using  the  polynomial  approach,  the  cy¬ 
clic  properties  can  be  utilized.  The  generator  polynomial  has  the  form 
g(x)  =  g0  +  g|X  +  g2x2  +  ...  +  gn.k  xn  k ,  and  a  block  of  information  data  can  be  represented 
as  a  polynomial  of  the  form  d(x)  =  d0  +  d,x  +  d2x2  +  ...  +  dk.,  xk~' .  Therefore,  a  codeword 
C  is  encoded  as  c(x)  =  d(x)g(x)  and  is  represented  as  c(x)  =  c0  +  c,x  +  c2x2  + ...  +  c^x"'1. 
Typically  RS  codes  are  put  into  a  systematic  format  in  order  to  reduce  the  number  of 
steps  needed  for  encoding.  A  codeword  is  in  systematic  format  if  it  has  the  form 
C  =  (p0,  p„  p2, ...  p2t_, ,  d0,  d„  d2, ...  dk., )  ,  where  the  data  symbols  (dj)  are  grouped  to¬ 
gether  and  the  n-k  parity  symbols  (pj)  are  also  grouped  together.  Using  the  cyclic  proper¬ 
ties  and  polynomial  representations,  the  codeword  has  the  form  c(x)  =  p(x)  +  x"'kd(x) , 
where  multiplying  by  x"'k  simply  performs  a  cyclic  shift  of  the  data  symbols  into  their 
proper  position.  In  the  above  polynomial  representation  the  parity  information  is 
equivalent  to  the  remainder  of  x"'kd(x)/g(x),  that  is  p(x)  =  x"‘kd(x)  mod  g(x)  [Has].  The 
resultant  code  word  polynomial  c(x),  systematic  or  non-systematic,  has  as  roots  the  same 
2t  consecutive  powers  of  a  as  does  the  generator  polynomial.  This  property  allows  for  an 
easy  way  of  determining  if  a  valid  codeword  has  been  received  and  for  efficient  decoding 
algorithms  of  RS  codes. 
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4.4  Decoding 


Since  the  introduction  of  RS  codes  in  1960,  there  has  been  a  rush  by  researchers  to 
find  an  efficient  decoding  algorithm  for  these  codes.  The  simple  techniques  used  in  lin¬ 
ear  block  codes  just  do  not  work.  For  instance  using  the  syndrome  look-up  table  method 
is  incomprehensible  due  to  the  number  of  syndrome  computations  needed.  For  instance, 
a  (63,53)  RS  code  requires  approximately  1020  syndrome  calculations  and  has  approxi¬ 
mately  1 095  (6453)  valid  codewords. 

In  the  original  paper  by  Reed  and  Solomon  [Ree]  they  introduced  a  possible  method 
for  decoding.  This  decoding  method  involved  the  solution  of  a  set  of  k  equations  in  k 
unknowns.  The  method  proved  to  only  be  useful  for  the  simplest  RS  codes.  Shortly  af¬ 
ter,  Peterson  [Pet2]  presented  a  popular  decoding  method  for  BCH  codes.  Gorenstein  and 
Zierler  [Gor],  and  later  Chien  [Chi]  and  Forney  [For2],  expanded  the  Peterson  algorithm 
for  BCH  codes  to  the  non-binary  case  for  RS  codes.  These  codes,  although  much  more 
efficient  than  any  others,  were  only  applicable  to  codes  with  a  small  number  of  errors. 

It  was  not  until  Berlekamp  [Ber2]  presented  his  version  of  an  efficient  decoder  for  use 
with  non-binary  BCH  and  RS  codes  that  the  application  of  RS  codes  became  more  than  a 
mathematical  exercise.  One  year  later  Massey  [Mas]  showed  that  the  Berlekamp  algo¬ 
rithm  could  be  realized  as  a  technique  of  finding  the  shortest  Linear  Feedback  Shift  Reg¬ 
ister  (LFSR)  capable  of  generating  a  given  sequence.  The  combination  of  the  two  ideas 
has  now  become  what  is  known  as  the  Berlekamp-Massey  algorithm.  Although  the  Ber- 
lekamp-Massey  algorithm  is  not  the  most  mathematically  simple  to  understand,  it  pro¬ 
duces  an  efficient  technique  for  decoding  RS  codes.  The  explanation  of  the  Berlekamp- 
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Massey  algorithm  that  follows  is  a  simple  overview  and  the  reader  is  referred  to  [Berl], 
[Ber3],  [Has],  [Lin],  [Nas],  [Wicl], 

The  Berlekamp-Massey  algorithm  is  a  four-step  process  for  error  correction  in  RS 
codes.  Reed-Solomon  codes,  which  unlike  BCH  codes  are  word  oriented,  require  an  ad¬ 
ditional  step  of  determining  the  error  value.  The  first  step  in  the  Berlekamp-Massey  algo¬ 
rithm  is  to  calculate  the  syndromes.  Next  the  error  locator  polynomial  must  be  obtained. 
Once  the  error  locator  polynomial  has  been  obtained,  the  error  location  can  be  found.  Fi¬ 
nally,  the  error  values  are  calculated.  These  four  steps  are  explained  in  the  following 
paragraphs. 

Recall  that  a  received  vector  has  the  form  r(x)  =  r0  +  r,x  +  r2x2  +  ...  +  rn.|Xn  I.  The  re¬ 
ceived  vector  can  be  represented  has  the  combination  of  a  valid  codeword  and  an  error 
pattern,  r(x)  =  c(x)  +  e(x)  where  c(x)  =  c0  +  c,x  +  c2x2  +  ...  +  cn.,xn''  is  a  valid  codeword 
and  e(x)  =  e0  +  e,x  +  e2x2  +  ...  +  en.,xn'1  is  an  error  pattern.  Since  the  codeword  c(x)  has  as 
roots  the  2t  consecutive  powers  of  a,  the  received  vector  will  have  the  same  roots  if  no 
errors  are  present.  That  is  r(x)  =  c(x)  when  e(x)  is  equal  to  all  zeros.  Calculating  the  syn¬ 
dromes  is  done  by  evaluating  the  received  vector  at  all  2t  consecutive  roots  of  the  code¬ 
words,  Sj  =  r(a')  =  c(a‘)  +  e(a')  for  i=  l,2,...,2t.  Since  S;  will  equal  zero  except  where 
e(a')  is  not  zero,  the  syndrome  is  calculated  from  the  equation  below  to  find  out  if  errors 
exist. 

n-1 

Si=H  ek(a')k>  i  =  l,2,...,2t 
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The  second  step  of  the  Berlekamp-Massey  algorithm  is  to  obtain  the  error  locator 
polynomial  A(x).  As  the  name  implies,  this  polynomial  gives  information  as  to  the  loca¬ 
tions  of  the  errors  in  the  received  code- vector.  It  is  actually  the  integer  power  of  the  in¬ 
verse  root  of  the  locator  polynomial  that  reveals  the  error  location.  Simply  stated:  find 
the  roots  of  the  error  locator  polynomial,  find  the  inverse  of  the  roots,  and  then  the  power 
of  the  resultant  is  the  location  of  the  error.  For  example  if  the  root  was  a2  defined  over 
GF(8),  then  the  inverse  of  a2  is  a5  and  the  error  is  at  the  fifth  location. 

It  is  the  method  of  finding  the  error  locator  polynomial  that  makes  the  Berlekamp- 
Massey  algorithm  so  complex.  The  error  locator  polynomial  is  found  through  an  iterative 
process,  usually  performed  in  tabular  form,  whereby  the  solution  for  A(x)  is  checked  for 
correctness  after  each  iteration.  If  the  error  locator  polynomial  A(x)  is  not  correct,  then  a 
correction  factor  is  added  to  A(x)  and  the  process  continues  for  2t  iterations.  The  error 
locator  polynomial  A(x)  is  given  below  for  v  errors,  where  A{  is  a  coefficient  of  the  form 
ot’. 

A(x)  =  (1+Aj  X+A2  X2+...+Av  Xv) 

After  the  error  locator  polynomial  A(x)  has  been  found,  the  third  step  is  to  find  all  of 
the  roots  of  the  locator  polynomial.  This  can  be  performed  by  factoring  A(x)  into  the 
form  A(x)  =  (l+AiXi)(l+Aj+1X'+1)...  (1+AVXV)  or  just  by  evaluating  A(x)  at  all  power  of 
a.  Once  the  roots  of  A(x)  are  found,  the  inverses  of  the  roots  reveal  the  error  locations. 

The  fourth  step  in  the  Berlekamp-Massey  algorithm  is  that  of  calculating  the  actual 
error  values  at  the  now  known  error  locations.  In  order  to  perform  this  task,  the  polyno- 
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mial  known  as  the  error  magnitude  polynomial  Z(x)  is  defined  using  the  error  locator 
polynomial  coefficients  and  the  syndrome  values.  The  error  magnitude  polynomial  is 
defined  as  Z(x)  =  1  +  (S,  +  A,)x  +  (S,  +  A,  S,  +  A2)  x2  +...+  (Sv  +  A,  Sv.,  +...+  Av)  xv  for 
v  errors.  This  error  magnitude  polynomial  is  used  in  the  following  equation  to  find  the 
error  values.  The  numerator  of  the  error  value  equation  is  the  error 

Z(X-*) 

ey =  — ; - 

n  o+XiX-V) 

i  =1 

'*j 

magnitude  polynomial  evaluated  at  the  root  of  the  error  locator  polynomial.  The  de¬ 
nominator  is  equal  to  v  factors  where  X,  represents  all  roots  of  the  error  locator  polyno¬ 
mial  that  are  not  at  the  current  location  being  calculated,  and  X'j  is  the  inverse  root  of  the 
error  locator  polynomial  for  the  current  error  value.  Once  all  v  values  for  e^  are  found  the 
error  polynomial  can  be  constructed  and  added  back  to  the  received  vector  for  error  cor¬ 
rection.  The  error  polynomial  would  have  the  form  e(x)  =  ej  +  ei+1x,+1  +  ei+2x,+2  +  ...  +  en. 
,xn'1  and  the  corrected  codeword  are  obtained  by  c(x)  =  r(x)  +  e(x). 

5  Interleaving 

Communication  systems  discussed  up  to  this  point  have  been  shown  to  operate  over 
channels  with  memory.  Furthermore,  it  is  the  memory  in  the  channel  that  is  to  blame  for 
the  majority  of  burst  errors.  Interleaving  is  a  technique  that  can  be  applied  in  an  attempt 
to  randomize  the  transmission  of  encoded  signals  over  the  bursty  channel.  This  random¬ 
izing  of  the  data  allows  the  channel  to  appear  memory  less.  The  remainder  of  this  section 
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will  discuss  the  two  most  widely  applied  interleaving  techniques,  block  interleaving  and 
convolutional  interleaving. 

5.1  Block  Interleaving 

Block  interleaving  attempts  to  scramble  the  codeword  over  the  channel  so  that  the 
channel  appears  to  be  random.  Block  interleaving,  although  appearing  random  to  the 
channel,  is  truly  a  one-to-one  deterministic  ordering  of  the  codeword  symbols  [Ram]. 

This  one-to-one  sequencing  of  the  data  provides  for  a  very  efficient  process  in  that  un¬ 
scrambling  of  codewords  is  easily  carried  out.  The  unscrambling,  or  deinterleaving  as  the 
process  will  be  called  hereafter,  is  essentially  an  interleaver  used  in  reverse  order.  All  of 
the  important  properties  will  be  brought  forth  in  the  following  discussion  of  the  block  in¬ 
terleaving  process. 

A  block  interleaver  is  nothing  more  than  a  buffering  mechanism  inserted  between  the 
error  coding  circuit  (channel  encoder)  and  the  modulator.  Similarly,  the  de-interleaver  is 
a  buffer  placed  after  the  demodulator  and  before  the  channel  decoder.  Figure  5.0  is  a 
simplified  block  diagram  of  a  typical  communication  system.  The  interleaver  and  de¬ 
interleaver  can  be  thought  of  as  a  d  x  n  matrix  where  d  is  the  depth  of  the  interleaver  indi¬ 
cating  the  number  of  codewords  stored.  Each  cell  in  the  matrix  represents  one  symbol  of 
length  m  bits.  Depending  on  the  size  of  the  interleaver  desired,  it  is  normally  imple¬ 
mented  with  memory  devices.  Therefore,  in  the  case  of  RS  codes,  memory  devices  of 
n-m-d  bits  are  required. 
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Figure  5.0  Simplified  Communications  System  Block  Diagram 


The  operation  of  the  block  interleaving  process,  as  illustrated  in  Figure  5.1,  is  simply 
the  inputting  of  codewords  into  the  interleaver  matrix  in  a  row-by-row  fashion,  then  out- 
putting  them  to  the  modulator  column-by-column.  This  technique  acts  to  spread  the 
codewords  over  the  channel  and  therefore  spreads  the  effects  of  a  burst  error  over  multi¬ 
ple  codewords.  With  block  interleaving,  as  shown  in  Figure  5.2,  a  burst  error  of  length  b 
<  d  in  the  binary  case  or  b  <  m  d  in  the  non-binary  case,  will  be  reduced  to  a  single  er¬ 
ror  in  each  codeword  [Ste].  Therefore,  the  length  of  the  burst  error  to  be  corrected  is  de¬ 
pendent  on  the  depth  (d)  of  the  block  interleaver.  The  depth  of  the  block  interleaver  can 
be  increased  but  not  without  the  tradeoff  of  increased  delays.  Also  as  shown  in  Figure 
5.3,  any  burst  error  of  length  b  =  r  m-d  where  r  >  1,  will  produce  r  symbol  errors  in  each 
codeword.  Hence,  a  code  with  error-correcting  ability  t  >  r  would  be  needed.  Therefore, 
by  using  block  interleaving  techniques,  a  less  complex  of  the  error  correction  code  can  be 
applied.  Reed-Solomon  codes  with  lower  error-correcting  capabilities  can  be  used  in 
conjunction  with  interleaving  techniques  to  correct  large  burst  errors. 
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Input  Stream:  A1 ,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  A1 0,  All ,  A1  2,  B1 ,  B2,  B3,  B4,  B5,  B6, . ,F9,  F10.  FI  1,  F12 


Output  Stream:  A1 3 1 ,  Cl ,  D1 ,  El ,  FI ,  A2,  B2,  C2,  D2,  E2,  F2,  A3,B3,  C3,  D3,  E3,  F3,  A4, . ,  C12,  D12,  El  2,  FI  2 

Figure  5.1  Block  Interleaver 

Conversely,  using  a  RS  code  with  a  higher  error  correction  capability  and  applying  inter¬ 
leaving,  further  strengthens  the  code  and  allows  for  the  correction  of  much  longer  burst. 
Note,  all  data  in  the  following  figures  are  read  from  left  to  right  with  the  left  most  data 
element  being  the  first  element  input  to  the  interleaver/de-interleaver.  The  shaded  rectan¬ 
gles  represent  the  occurrence  of  symbol  errors. 
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Interleaver 


Input  Stream:  A1 ,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  A 10,  A1 1 ,  A 1 2,  B1 ,  B2,  B3,  B4,  B5,  B6, 


,,F9,F10,  FI  1 ,  F12 
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Output  Stream: 


E2,  F2,  A3,  B3,  C3,  D3,  E3,  F3,  A4, . C12,  D12,  E12,  F12 


burst 


De-interleaver 


Input  Stream:  A1,B1,  Cl,  Dl,  g|§ 


E2,  F2,  A3,  B3,  C3,  D3,  E3,  F3,  A4, . 
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Output  Stream:  A1 


A1 

f  A2 

A3,  A4,  A5,  A6,  A7,  A8,  A9,  A1 0,  A1 1 ,  A1 2,  B1 , 

ii 

§§f 

i  ii 

B3,  B4,  B5,  B6, . JF9,F10,  FI  1 ,  F12 


Figure  5.2  Block  Interleaver  with  (b  <  d) 
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Interleaver 


Input  Stream:  Al,  A2,  A3,  A4,  A5,  A6,  A7,  A8,  A9,  A10,  A1 1,  A12,  Bl,  B2,  B3,  B4,  B5,  B6, . ,F9,F10,  FI  1,  F12 
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Output  Stream:  r 


A3,  B3,  C3,  D3,  E3,  F3,  A4, . ,C12,  D12,  E12,  F12 


burst 


De-interleaver 


Input  Stream: 


A3,B3,C3,D3,E3,F3,  A4, 


4 _  n  symbols  - ► 


B 

■ 

a 

a 

g 

1^3 

p 

P 

a 

PI 

ggiii 

j|§ 

H 

yi 

■ra 

gmi 

(g 

H 

p 

a 

|| 

m 

a 

p| 

pp 

m 

H 

a 

a 

a 

a 

p 

n 

B 

gg 

m 

g 

a 

a 

g 

p 

p 

H 

B 

p 

a 

a 

a 

p 

g 
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Output  Stream:  fAl,A2l  A3,  A4,  A5,  A6,  A7,  A8,  A9,  A10,  All,  A12,  B3,  B4,  B5,  B6, 


,,F9,F10,  F11.F12 


Figure  5.3  Block  Interleaver  with  (b  =  rd;  r  =  2) 
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5.2  Convolutional  Interleaving 


Convolutional  interleavers  are  very  much  like  block  interleavers  in  that  their  function 
is  to  scramble  the  codewords  before  transmission.  The  scrambling  of  a  convolutional  in¬ 
terleaver  is  also  a  one-to-one  deterministic  process.  It  can  be  seen  that  for  a  block  inter¬ 
leaver,  after  a  column  of  d  code  symbols  has  been  output  from  the  interleaver,  storage 
space  is  still  allocated  for  the  already  transmitted  column  of  data.  The  additional  storage 
needed  for  the  block  interleaver  produces  unwanted  delays  and  memory  requirements  in 
the  communication  system.  Convolutional  interleaving  is  a  more  optimal  interleaving 
technique  in  that  it  reduces  the  amount  of  storage  space  required  and  therefore  reduces 
delays.  There  are  many  excellent  references  available  describing  interleaving  techniques, 
however  the  explanation  of  convolutional  interleaving  that  follows  is  as  proposed  by 
[Fori]  and  further  explained  by  [Ste]  and  [Siml]. 

Figure  5.4  is  a  representation  of  a  convolutional  interleaver  and  corresponding  de¬ 
interleaver.  The  first  input  line  of  the  interleaver  is  connected  directly  to  the  channel 
(modulation  is  assumed).  All  other  input  lines  are  composed  of  registers  with  each  suc¬ 
cessive  register  containing  £  additional  storage  cells.  A  serial  data  stream  of  codewords 
from  the  encoder  circuit  is  input  to  the  convolutional  interleaver.  Serial  to  parallel  con¬ 
version  maybe  needed.  The  second  input  line  is  input  to  a  shift  register  of  length  £  sym¬ 
bols.  The  third  line  is  input  into  a  shift  register  of  length  2  £  symbols  and  the  process 
continues  for  d-1  lines,  for  a  total  of  d  input  lines.  Remember  that  for  RS  codes  each 
symbol  is  m  bits  in  length  and  therefore,  each  shift  register  cell  is  considered  to  be  m  bits 
in  length.  At  the  same  time  that  the  input  demultiplexer  is  inputting  a  codeword  symbol. 
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the  output  multiplexer  is  selecting  the  oldest  symbol  shifted  out  from  each  shift  register. 
The  output  stream  is  then  sent  to  the  modulator. 


Figure  5.4  Convolutional  Interleaver 


For  the  case  of  using  convolutional  interleavers  to  scramble  RS  coded  information,  it  is 
convenient  to  design  the  interleaver  with  the  depth  (d)  equal  to  the  length  of  a  codeword 
(d  =  n).  Therefore,  a  delay  of  d  l  symbols  or  d  im  bits  is  placed  between  adjacent  input 
symbols.  From  this  implementation,  it  is  obvious  that  a  burst  of  length  (d  ^+l)m  bits  will 
only  introduce  one  symbol  error  in  d  t  consecutive  codewords  after  de-interleaving  takes 
place. 


5.3  Pseudorandom  Interleaving 

Another  important  interleaving  technique  worth  mentioning  is  called  pseudorandom 
interleaving.  Pseudorandom  interleavers  work  on  the  principle  of  using  synchronous 
pseudorandom  number  generators  in  both  the  interleaver  and  deinterleaver.  The  pseu¬ 
dorandom  codes  are  used  to  vary  the  depth  of  interleaving  in  a  manner  only  known  to  the 
communication  system.  The  importance  of  this  technique  is  in  combating  intentional 
jamming  in  hostile  communication  environments.  If  the  depth  of  an  interleaver  is  known. 
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a  pulsed  jammer  can  be  used  with  a  pulse  repetition  interval  (PRI)  equal  to  the  interleav¬ 
ing  depth  in  order  to  produce  burst  errors  at  the  output  of  the  deinterleaver  [Sim2].  By 
using  pseudorandom  interleaving  depths,  pulsed  jammers  lose  their  effectiveness.  The 
information  provided  here  on  pseudorandom  interleaving  is  only  for  completeness  and 
will  not  be  explored  any  further. 

5.4  Interleaving  Delay 

Referring  back  to  Figure  5.1,  the  block  interleaving  scheme  uses  d  •  n  symbol  storage 
cells  or  d  •  m  •  n  bit  storage  cells.  This  introduces  a  delay  of  n(d-l)+l  symbols  because 
transmission  can  begin  as  soon  as  the  first  symbol  in  the  last  row  is  input.  Obviously, 
this  allows  for  the  output  of  the  first  column  before  the  input  of  the  remaining  symbols  in 
the  last  row.  If  measuring  the  delay  of  the  interleaver  in  bits,  then  there  is  a  m[n(d-l)+l] 
bit  delay.  The  delay  can  be  measured  in  time  (seconds)  using  the  known  fact  that  the  data 
rate  (  R )  is  equal  to  the  inverse  of  the  bit  transmission  rate  (R  =  l/tg).  Therefore,  the  time 
delay  is  (ts*TB)  seconds  where  TB  =  m[n(d-l)+l]  is  the  interleaver  bit  delay  measured  by 
the  number  of  bits.  Therefore,  the  end-to-end  delay  time  for  block  interleaving  is 
2(ts‘TB)  +  channel  delays. 

Convolutional  interleavers  are  more  efficient  than  block  interleavers  because  of  their 
smaller  storage  requirement  and  therefore  smaller  delay  times.  Referring  to  Figure  5.4,  it 
can  be  seen  that  the  total  number  of  storage  cells  needed  for  a  convolutional  interleaver  is 
n-  £( d- 1  )/2  symbols.  For  simplicity  define  £  to  be  equal  to  one  and  recall  each  symbol  is 
m  bits  in  length.  Therefore,  the  storage  requirement  can  be  measured  as  m[n(d-l)]/2  bits. 
As  it  was  done  for  the  block  interleaving,  time  delay  can  be  measured  as  ts-m*[n(d-l)]/2 
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seconds.  Recalling  that  for  the  block  interleaver  the  storage  requirement  is  n(d-l)+l 
symbols  and  n(d-l)/2  symbols  for  the  convolutional  interleaver,  therefore  it  can  be  seen 
that  the  convolutional  interleaver  is  nearly  twice  as  efficient  as  the  block  interleaver.  The 
interleaver  efficiency  is  examined  in  more  detail  in  [Wicl], 
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6  Interleaved  Burst  Error  Correction  Efficiency 

Up  to  this  point  a  preliminary  review  of  RS  codes  and  error  control  coding  in  general 
have  been  presented.  This  section,  along  with  Section  7,  constitutes  the  findings  of  the 
research  performed  on  evaluating  the  burst  error  correction  of  RS  codes.  First  in  this 
section  the  burst  error  correction  efficiency  of  linear  block  codes  will  be  explained  using 
the  traditional  definition  of  efficiency.  The  concept  of  phased  burst  error  correction  effi¬ 
ciency  will  also  be  explained.  Then  finally  the  burst  error  correction  efficiency  of  RS 
codes  will  be  defined. 

Recall  from  Section  2  one  that  burst  errors  are  defined  according  to  the  length  of  the 
burst  b,  and  the  guard  space  g  separating  the  individual  burst  errors.  Hence,  the  Gallager 
bound  states  that  for  a  binary  code  of  burst  error-correcting  capability  band  guard  space 
g,  the  code  rate  R  must  satisfy  the  inequality  below  [Gal2].  The  Reiger  bound,  which  is 


derived  from  the  Gallager  bound,  states  that  a  (n,k)  linear  block  code  with  burst  error- 
correcting  capability  b  has  to  meet  the  following  condition  [Cap]. 

n  -  k  >  2b 

Also,  a  linear  block  code  with  burst  error  correction  b  and  detection  d  must  meet  the  con¬ 
dition  below. 

n  -  k  +  \  >  b  +  d 
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Based  on  the  Reiger  bound,  the  traditional  method  for  defining  burst  error  correction  effi¬ 
ciency  is  z  as  shown  below.  This  ratio  is  called  the  “phased  burst  error  correction  effi¬ 
ciency”.  Linear  block  codes  for  which  z  =  1  are  considered  optimal  codes  because 


these  give  the  most  error  correction  achievable  for  the  amount  of  parity  used. 

6. 1  Maximum  Burst  Error  Correction  of  Reed-Solomon  Codes 

The  Gallager  bound,  Reiger  Bound,  and  the  phased  burst  error  correction  efficiency 
are  defined  for  binary  linear  block  codes  correcting  single  burst  errors.  Reed-Solomon 
codes  are  non-binary  multiple  burst  error-correcting  linear  block  codes.  Recall  from  Sec¬ 
tion  4  that  RS  codes  have  random  error-correcting  capability  t  =  (n-k)/2  and  maximum 
single  burst  error-correcting  capability  bImax=  tm,  where  the  subscript  (1)  indicates  single 
burst  error  correction.  The  term  maximum  is  used  because  the  burst  errors  must  be  lo¬ 
cated  exactly  over  t  symbols  (which  may  not  always  be  the  case).  This  (blmax)  forms  an 
upper  bound  on  the  single  burst  error-correcting  capability  of  RS  codes.  Using  the  phased 
single  burst  error-correcting  capability  of  RS  codes  blmax  and  t  =  (n-k)/2,  it  can  be  seen 
that  RS  codes  are  always  optimal  phased  single  burst  error-correcting  codes. 


zi 


2  b ,  2  tm 

(n  -  k )  ( n  -  k  )  m 


(n  -  k  )  m 


( n  -  k  )  m 
(n  -  k ) m 
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6.2  Multiple  Burst  Error  Correction  of  Reed-Solomon  Codes 


In  some  situations  frequent  bursts  occur  in  a  given  codeword  and  single  burst  error- 
correcting  codes  are  of  no  use.  This  usually  occurs  when  the  guard  spaces  are  not  well 
defined  and  when  a  cluster  of  bursts  are  too  long  to  be  treated  as  a  single  large  burst.  As 
mentioned  earlier  RS  codes  are  multiple  burst  error-correcting  codes,  which  has  made 
them  very  popular  in  combating  these  types  of  errors.  With  this  in  mind,  the  phased 
(maximum)  double  and  triple  burst  error-correcting  capability  of  RS  codes  are  given  be¬ 
low.  Where  the  function  L*J  is  the  floor  function  and  is  defined  as  the  largest  integer  less 
than  or  equal  to  x.  Because  RS  codes  are  word-oriented  codes,  the  term  phased  burst 


error  correction  applies  only  to  the  case  that  the  burst  errors  occur  over  exactly  t  symbols. 
Therefore,  the  phased  burst  error  correction  of  RS  codes  is  the  maximum  burst  error  cor¬ 
rection. 

The  maximum  double-burst  error  correction  capability  is  defined  as  the  maximum  size 
of  two  simultaneous  burst  errors  that  can  be  corrected.  The  code  may  be  able  to  correct 
different  combinations  of  double  burst  errors  so  long  as  t  is  not  exceeded.  However,  no 
two  burst  greater  than  b2max  can  be  corrected  at  once.  For  instance  if  m  =  8  and  t  =  6 
((255,243)  RS  code),  the  code  will  correct  two  burst  of  maximum  length  24  bits.  How¬ 
ever,  the  same  code  may  correct  two  bursts,  one  of  length  16  bits  and  the  other  32  bits.  If 
one  burst  is  of  length  greater  than  b2max,  then  the  other  burst  will  be  less  than  b2max.  The 
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maximum  burst-error  correction  capability  also  applies  to  triple  and  higher  order  burst 
error-correcting  codes.  The  general  form  for  phased  burst  error  correction  efficiency  of  a 
RS  codes  is  given  as  zs  below.  It  can  be  easily  seen  that  although  the  double,  triple,  and 
higher  order  phased  burst  error  correction  efficiencies  of  RS  codes  are  not  optimal,  they 
have  very  high  efficiencies. 


where; 


z„  = 


min  -  k ) 


6.3  Guaranteed  Burst  Error  Correction  Of  Reed-Solomon  Codes 

The  maximum  or  phased  burst  error  correction  capabilities  defined  above  are  only  ap¬ 
plicable  when  the  burst  errors  occur  exactly  over  t  symbols.  In  most  practical  applica¬ 
tions  the  burst  errors  will  not  occur  exactly  over  t  symbols  and  therefore  the  need  to  find 
the  guaranteed  burst  error-correcting  capabilities  of  a  code.  Recall  from  Section  4  that 
the  guaranteed  single,  double,  and  triple  burst  error  correction  capability  of  RS  codes  are 
b„  b2,  and  b3  restated  below.  This  is  the  burst  error-correcting  capability  of  a  RS  code  no 
matter  where  the  burst  errors  occur. 

b ,  =  m  ( t  -  1 )  +  1 


299 


At  this  point  the  guaranteed  burst  error-correcting  efficiency  for  a  given  RS  code  is 
defined.  Using  the  traditional  definition  of  efficiency  and  the  guaranteed  burst  error- 


correcting  capability  equations,  the  single,  double,  and  triple  guaranteed  burst  error- 
correcting  efficiency  of  RS  codes  are  y„  y2  andy3  as  follows. 

2[w(/-l)  +  l] 
m(n  -  k) 
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-1 
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+  1 

m(n  -  k) 


r  f 
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-1 
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+  1 

m(n  -  k  ) 


For  the  general  case,  the  burst  error  correction  efficiency  of  an  s  burst  error-correcting  RS 
code  is  given  [Kam3]. 


where; 


2  sbs 

"^s  m{n-k) 


+  1 


6.4  Analysis 

Reed-Solomon  codes  of  length  n  =  255,  symbol  length  m  =  8  bits,  are  very  popular 
and  find  many  applications  due  to  their  byte  sized  symbols.  Also,  codes  with  rate  R  less 
than  fifty  percent  are  not  often  used  because  of  the  large  amount  of  redundancy.  There¬ 
fore,  Table  6.0  shows  RS  single,  double,  and  triple  guaranteed  burst  error  correction  ca- 
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pabilities  and  efficiencies  for  codes  of  length  n  =  255  and  rate  R  >  50%.  For  a  complete 
table  of  RS  codes  of  length  n  =  7, 15, 31, 63,  127,  and  255  see  appendix  A. 

As  can  be  seen  from  Table  6.0  and  as  illustrated  in  Figure  6.0,  the  guaranteed  single 
burst  error  correction  efficiency  approaches  optimality  very  quickly  as  code  rate  de¬ 
creases.  Note  that  as  can  be  observed  from  the  equations  for  guaranteed  burst  error  cor¬ 
rection  efficiency  of  RS  codes,  y, ,  y2 ,  and  y3  can  never  equal  one,  however  they  are  re¬ 
markably  close.  The  fluctuations  in  y2  and  y3  are  due  to  the  fact  that  there  are  RS  codes 
that  have  the  same  double  and  triple  guaranteed  burst  error-correcting  capabilities  with 
increased  redundancy  as  observed  in  the  flat  parts  in  the  graphs  of  Figure  6.1.  Figure  6.1 
only  shows  RS  codes  up  to  a  code  rate  of  eighty  percent  for  better  resolution.  For  a  com¬ 
plete  set  of  graphs  of  burst  error-correcting  capabilities  of  RS  codes  see  appendix  A. 

Using  the  guaranteed  burst  error-correcting  capability  and  efficiency  of  RS  codes,  a 
code  designer  now  has  additional  information  at  his/her  disposal  for  selecting  the  best 
possible  RS  code  for  the  error  correction  scheme  being  designed.  For  instance,  if  the  de¬ 
signer  is  interested  in  double  burst  error-correcting  codes,  the  codes  on  the  flat  areas  of 
the  guaranteed  double  burst  error-correcting  capability  graph  with  the  highest  code  rate 
are  the  best  choice.  The  other  codes  on  each  flat  region  have  decreased  code  rates  and 
lower  burst  error-correcting  efficiencies.  The  information  presented  here  can  help  guide 
the  code  designer  to  get  more  error  correction  for  the  amount  of  redundancy  (more  bang 
for  the  buck). 
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Table  6.0  Burst  Error  Correction  of  Reed-Solomon  Codes  of  Length  n  =  255 
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Burst  Error  Correcting  Capabiiity  (bits) 


BURST  ERROR  CORRECTION  EFFICIENCY 
(Reed-Solomon  Codes  of  Length  n  «  255) 


-Single  Burst 
-Double  Burst 
-Triple  Burst 


Figure  6.0  Guaranteed  Burst  Error  Correction  Efficiency 


GUARANTEED  BURST  ERROR  CORRECTING  CAPABILITY 
(Reed-Solomon  Codes  of  Length  n  =  255) 
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Figure  6.1  Guaranteed  Burst  Error  Correction  Capability 
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7  Burst  Error  Correction  Efficiency 


Up  to  this  point  the  classical  definition  of  efficiency  for  a  linear  block  code  has  been 
expanded  to  include  non-binary  RS  codes.  This  research  has  produced  efficiency  equa¬ 
tions  that  will  be  useful  in  designing  RS  coding  schemes.  However,  due  to  the  fact  that 
interleaving  is  commonly  used  in  combating  bursty  errors,  a  look  at  the  efficiency  of  in¬ 
terleaved  RS  codes  is  in  order.  The  remainder  of  this  thesis  will  investigate  the  efficiency 
of  interleaved  RS  codes.  The  question  of  whether  or  not  the  efficiency  of  a  RS  coding 
scheme  increases  with  the  use  of  interleaving  will  be  answered. 

7.1  Block  Interleaved  Guaranteed  Burst  Error  Correction 

Recall  that  block  interleaving  scrambles  codewords  to  be  transmitted  over  the  com¬ 
munications  channel  thus  randomizing  the  effects  of  a  burst  error.  Essentially,  the  burst 
error  is  broken  into  smaller  and  more  manageable  errors  over  multiple  codewords.  It  is 
the  depth  (d)  of  the  block  interleaver  that  determines  how  burst  errors  will  be  randomized 
over  the  interleaved  codewords.  Since  the  depth  of  the  block  interleaver  is  equal  to  the 
number  of  codewords  being  interleaved,  a  burst  of  length  b  =  rdm  bits  (where  r  =  2,3,4,... 
is  the  number  of  columns  in  the  interleaver  with  errors,  see  Figures  5.2  and  5.3)  will  pro¬ 
duce  burst  errors  of  b  =  rm  bits  maximum  in  each  de-interleaved  codeword.  Therefore, 
codewords  with  the  ability  to  correct  burst  errors  of  length  b  =  rm  bits  and  an  interleaver 
of  depth  d  is  all  that  is  required  to  correct  burst  of  length  b  -  rdm  bits. 

Reed-Solomon  codes  that  can  correct  t  =  r  random  symbols  errors  can  be  used  to  cor¬ 
rect  burst  errors  up  to  b  =  rm  bits.  However,  there  are  instances  when  burst  errors  do  not 
occur  exactly  over  t  =  r  symbols  and  the  random  error-correcting  capability  t  is  exceeded 
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(t  <  r).  The  equations  for  guaranteed  burst  error  correction  were  given  in  Section  4  and 
are  restated  here.  Guaranteed  single  burst  error  correction  is  b;  =  m(t  - 1)  + 1  bits,  guar¬ 
anteed  double  burst  error  correction  is  b2  =  m(j_t/2j  - 1)  +  1  bits  and  guaranteed  triple 
burst  error  correction  is  b3  =  m(Lt/3j  - 1)  +  1  bits.  For  the  general  case  of  s  burst  errors 
the  equation  is  bs  =  m(Lt/sJ  - 1)  +  1.  These  equations  again  give  the  burst  error  correction 
of  a  given  RS  code  no  matter  where  the  burst  error  occurs  in  the  codeword. 

Combining  the  operation  of  the  block  interleaver  with  the  guaranteed  burst  error- 
correcting  capability  of  a  given  RS  code,  the  following  equations  are  derived. 


y/x  =  m(td  - 1)  + 1 


m 

t 

d- 1 

.2. 

+  1 


Vi  =m 


d- 1 


+  1 


\d- 1 


+  1 


LLSJ  J 

These  equations  express  the  single,  double,  and  triple  guaranteed  burst  error-correcting 
capability  of  a  block  interleaved  Reed-Solomon  code,  where  \|/s  is  the  general  form  of  the 
equation.  These  equations  simply  mean  that  using  a  known  RS  code  (n,k,t,m)  and  a 
block  interleaver  of  known  depth  (d),  a  burst  error  of  length  \j/  is  guaranteed  to  be  cor¬ 
rected  no  matter  where  it  occurs  over  the  interleaved  codewords. 
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7.2  Block  Interleaved  Guaranteed  Burst  Error  Correction  Efficiency 


In  Section  6  it  was  stated  that  from  the  Reiger  bound  and  the  Gallager  bound  ,  the 
phased  burst  error  correction  efficiency  of  a  linear  block  code  is  defined  as  z  =  2b/(n-k). 

If  z  =  1 ,  then  the  linear  block  code  is  said  to  be  optimal  giving  the  most  burst  error  cor¬ 
rection  for  the  amount  of  parity.  This  definition  of  efficiency  was  for  binary  linear  block 
codes  which  was  expanded  for  non-binary  RS  codes  later  in  Section  6  of  this  thesis 
(see  y, ,  y2,  y3 ,  and  ys ,  in  Section  6). 

Taking  the  definition  of  efficiency  (z)  another  step  forward,  it  will  now  be  expanded  to 
give  the  efficiency  of  block  interleaved  RS  codes.  Since  the  block  interleaver  increases 
the  number  of  parity  bits  being  transmitted  to  achieve  the  guaranteed  burst  error  correc¬ 
tion  given  in  \j/,  the  depth  (d)  is  applied  to  the  denominator  of  z.  As  before  m  is  applied 
to  the  denominator  because  RS  codes  are  non-binary.  Replacing  the  burst  error  correc¬ 
tion  factor  b  in  the  numerator  of  z  with  y,  the  equations  for  the  block  interleaved  guaran¬ 
teed  burst  error  correction  efficiency  of  a  RS  code  have  been  derived  and  are  listed  below. 
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The  importance  of  the  above  efficiency  equations  can  be  summarized  best  by  exam¬ 
ining  Figure  7.0  and  Table  7.0  below.  Figure  7.0  is  a  plot  of  the  block  interleaved  guar¬ 
anteed  single  burst-error  correction  efficiency  ( vp, )  of  length  255  for  RS  codes  of  vary¬ 
ing  interleaving  depths.  Length  255  is  chosen  here  because  of  the  popular  use  of  the 
eight  bit  (m  =  8 )  symbols.  Also,  only  the  high  code  rate  (R  >  75%)  RS  codes  are  shown 
to  reduce  the  congestion  of  the  graph. 

It  is  immediately  obvious  from  Figure  7.0  that  an  increase  in  interleaver  depth  pro¬ 
duces  an  increase  in  efficiency  of  any  given  RS  code.  Even  more  importantly,  examining 
the  data  from  Table  7.0,  it  is  seen  that  a  lower  code  rate  ( R  )  non-interleaved  code  can  be 
replaced  (based  on  efficiency)  with  a  higher  code  rate  ( R  )  interleaved  RS  code.  This 
means  that  a  simpler  RS  code  can  be  interleaved  and  produce  the  same  burst  error  correc¬ 
tion  and  burst  error  correction  efficiency  as  that  of  a  more  powerful  non-interleaved  RS 
code.  By  using  a  simpler  RS  code  there  is  a  reduction  in  hardware  cost  and  the  number 
of  computations  needed  for  encoding  and  decoding  is  reduced.  As  a  reminder,  interleav¬ 
ing  will  produce  delays  in  communication  and/or  data  recording  systems.  Therefore,  re¬ 
placing  a  more  complex  RS  code  with  a  simpler  interleaved  code  can  only  be  done  if  the 
system  can  tolerate  the  interleaving  delays  introduced. 
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Table  7.0  Block  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency 
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BLOCK  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255) 


Reed-Solomon  Code  (n,k) 


Figure  7.0  Block  Interleaved  Guaranteed  Single  Burst  Error  Correction  Efficiency 

7.3  Convolutional  Interleaved  Guaranteed  Burst  Error  Correction 

Convolutional  interleaving  is  very  much  like  block  interleaving  in  that  the  idea  is  to 
scramble  codewords  and  make  them  appear  random  to  the  communications  channel. 
Convolutional  interleaving  was  shown  in  Section  5  to  be  more  efficient  than  block  inter¬ 
leaving,  in  that  it  requires  less  storage.  For  the  convolutional  interleaver  it  is  the  depth 

(d)  and  the  number  of  storage  cells  ( l )  on  each  level  that  determine  how  burst  errors  will 
be  randomized.  From  Figure  5.4  it  can  be  seen  that  each  symbol  of  a  burst  error  will  be 
separated  by  d  £  symbols  after  convolutional  interleaving. 
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If  X  =  (x0,  x„  x2,  x3,  x4, ...)  is  a  codeword  applied  to  the  convolutional  interleaver  of 
Figure  5.4,  then  the  data  stream  to  be  transmitted  over  the  communications  channel 
would  be  that  of  Figure  7. 1 .  From  Figure  7. 1  it  is  clear  that  a  single  burst  of  length 
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Figure  7.1  Output  of  a  Convolutional  Interleaver 


b  =  m(df  +1)  bits  will  produce  a  single  random  error  in  each  RS  codeword  after  deinter¬ 
leaving.  Combining  this  fact  about  convolutional  interleaving  with  the  concept  of  guar¬ 
anteed  burst  error  correction  (see  guaranteed  burst  error  correction  capability  in  Section 
6)  the  following  equations  are  derived.  These  equations  express  the  single,  double,  and 
triple  guaranteed  burst  error  correction  capability  of  a  convolutional  interleaved  RS  code 
no  matter  where  the  burst  error  occurs,  where  <|>s  is  the  general  form  of  the  equation  for  s 
burst  errors. 
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7.4  Convolutional  Interleaved  Guaranteed  Burst  Error  Correction  Efficiency 

As  stated  in  Section  6  the  phased  burst  error  correction  efficiency  is  classically  defined 

as  z  -  2b/(n-k).  Throughout  this  thesis  the  classical  definition  of  efficiency  has  been 
modified  for  various  applications  of  RS  codes.  In  Section  6  the  efficiency  (z)  was  ex¬ 
panded  for  the  case  of  guaranteed  burst  error  correction  for  RS  codes  (y„  Y2,Y3,Ys>)-  In 
this  section  the  definition  of  efficiency  (z)  has  been  expanded  to  include  the  case  of  block 
interleaved  guaranteed  burst  error  correction  efficiency  of  RS  codes.  Similarly,  in  this 
section  the  definition  of  efficiency  (z)  is  now  expanded  to  include  the  case  of  convolu¬ 
tional  interleaved  guaranteed  burst  error  correction  efficiency.  The  equations  provided 
below  are  derived  by  substituting  <(>„  <j)2s  <J>3,  and  <j>s  for  b  in  the  equation  z  =  2b/(n-k),  m  is 
applied  to  the  denominator  because  RS  codes  are  non-binary.  The  term  d  i  +  1  is  applied 
to  the  denominator  because  there  must  be  d  t  +  1  codewords  processed  before  achieving 
the  guaranteed  burst  error  correction  of  <(>,,  <j>2,  <t>3,  and  <|>s .  These  equations  represent  the 
guaranteed  burst  error  correction  efficiency  achievable  using  convolutional  interleaving. 
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As  can  be  observed  in  Figures  7.1  through  7.4  and  Tables  7.1  through  7.4  the  effi¬ 
ciency  of  any  given  RS  code  increases  with  an  increase  in  convolutional  interleaving 

depth.  Also,  it  is  clear  that  as  the  additional  delay  per  each  set  of  registers  (£)  is  increased 
the  single  burst  error-correcting  efficiency  is  increased.  Most  importantly,  notice  that  as 
the  amount  of  convolutional  interleaving  (delay)  is  increased  (here  as  £  increases  in  fig¬ 
ures  7.1  to  7.4)  the  efficiency  is  very  high  for  the  simplest  RS  codes.  For  example  the 
(255,253)  RS  code  with  convolutional  interleaving  of  d  =  2  and  £  =  4  (figure  7.4)  has  a 
guaranteed  single  burst  error  correction  efficiency  of  more  than  ninety  percent  and  can 
correct  a  single  burst  of  sixty-five  bits  with  a  code  rate  (R)  greater  than  ninety-nine  per¬ 
cent.  Very  impressive! 
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Table  7.1  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (*=1) 
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Burst  Error  Correction  Efficiency 


CONVOLimONAL  INTERLEAVED  Rffi^SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Connection  Efficiency 
(Code  Length  n= 255,  1  =  1) 


(255,253)  (255,247)  (255,241)  (255,235)  (256,229)  (255,223)  (255,217)  (255,211)  (255,205)  (255,199) 


Reed-Sdomon  Code  (n,k) 


Figure  7. 1  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (£=1) 


Table  7.2  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 
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Burst  Error  Correction  Efficiency 


CONVOLUTIONAL  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255,  1  =  2) 


Reed-Sdomon  Code  (n,k) 


Figure  7.2  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (7=2) 
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Table  7.3  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (£=3) 
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CONVOLUTIONAL  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255,  1  =  3) 


Reed-Sol  omon  Code  (n,k) 


Figure  7.3  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (£=3) 


Table  7.4  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Eff.  (T=4) 
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Burst  Error  Correction  Efficiency 


CONVOLUTIONAL  INTERLEAVED  REED-SOLOMON  CODE 
Guaranteed  Single  Burst  Error  Correction  Efficiency 
(Code  Length  n  =  255,  1=4) 


Reed-Solomon  Code  (n,k) 


Figure  7.4  Convolutional  Interleaved  Guaranteed  Single  Burst  Error  Correction 

Efficiency  (7=4) 
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8  Conclusion 


Burst  errors  are  commonplace  in  today’s  digital  communication  and  recording  sys¬ 
tems.  Reed-Solomon  codes  are  known  to  be  the  best  at  combating  burst  errors  because 
they  are  non-binary  Maximum  Distance  Separable  (MDS)  codes.  It  was  the  intent  of  this 
research  effort  to  explore  the  burst  error-correcting  capabilities  and  efficiencies  of  Reed- 
Solomon  codes. 

This  thesis  discussed  the  phased  burst  error-correcting  capabilities  of  RS  codes.  The 
concept  of  guaranteed  burst  error  correction  was  then  examined.  The  idea  of  guaranteed 
burst  error  correction  was  then  expanded  to  include  the  popular  block  and  convolutional 
interleaving  techniques  illustrating  the  increase  in  guaranteed  burst  error  correction 
achievable  using  interleaving.  The  classical  definition  of  error  correction  efficiency  was 
then  expanded  to  include  the  guaranteed  burst  error  correction  of  RS  codes  and  inter¬ 
leaved  RS  codes. 

It  was  shown  that  RS  codes  are  very  efficient  at  correcting  burst  errors.  In  fact  many 
RS  codes  approach  optimal  efficiencies  while  correcting  considerably  large  burst  errors. 
All  of  the  equations,  tables,  and  graphs  generated  in  this  thesis  can  be  used  by  code  de¬ 
signers  as  additional  criteria  in  choosing  the  most  efficient  burst  error-correcting  system 
for  their  needs. 
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BURST  ERROR  CORRECTION  (bits) 


9  Appendix  A 
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BURST  ERROR  CORRECTION  (bits)  EFFICIENCY 
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(Reed-Solomon  Codes  of  Length  n  =  63) 
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CODE  RATE 


BURST  ERROR  CORRECTING  CAPABILITY 
(Reed-Solomon  Codes  of  Length  n  =  127) 


t 


Y 


(n,k) 

(255,253)  1  0.992157  ~ 

(255,251)  2  0.984314  9 

(255,249)  3  0,976471  17 

(255,247)  4  0.968627  25 

(255,245)  5  0.960784  33 

(255,243)  6  0.952941  41 

(255,241)  7  0,945098  49 

(255,239)  8  0.937255  57 

(255,237)  9  0.929412  65 

(255,235)  10  0.921569  73 

(255,233)  11  0,913725  81 

(255,231)  12  0.905882  89 

(255,229)  13  0.898039  97 

(255,227)  14  0.890196  105 

(255,225)  15  0.882353  113 

(255,223)  16  0.874510  121 

(255,221)  17  0.866667  129 

(255,219)  18  0.858824  137 

(255,217)  19  0.850980  14? 

(255,215)  20  0.843137  153 

(255,213)  21  0.835294  161 

(255,211)  22  0.827451  169 

(255,209)  23  0.819608  177 

(255,207)  24  0.811765  185 

(255,205)  25  0.803922  193 

(255,203)  26  0.796078  201 

(255,201)  27  0.788235  209 

(255,199)  28  0.780392  217 

(255,197)  29  0.772549  225 

(255,195)  30  0.764706  233 

(255,193)  31  0.756863  2aT 

(255,191)  32  0.749020  249 

(255,189)  33  0.741176  257 

(255,187)  34  0.733333  265 

(255,185)  35  0.725490  273 

(255,183)  36  0.717647  281 

(255,181)  37  0.709804  289 

(255,179)  38  0,701961  297 

(255,177)  39  0.694118  305 

(255,175)  40  0.686275  313 

(255,173)  41  0,678431  321 

(255,171)  42  0.670588  329 

(255,169)  43  0.662745  337 

(255,167)  44  0.654902  345 

(255,165)  45  0.647059  353 

(255,163)  46  0.639216  361 

(255,161)  47  0.631373  369 

(255,159)  48  0.623529  377 

(255,157)  49  0.615686  385 

(255,155)  50  0.607843  393 

(255,153)  51  0.600000  4?l" 

(255,151)  52  0.592157  409 

(255,149)  53  0.584314  4l7 

(255,147)  54  0.576471  425 

(255,145)  55  0.568627  433 

(255,143)  56  0.560784  441 

(255,141)  57  0,552941  449 

(255,139)  58  0.545098  457 

(255,137)  59  0.537255  465 

(255,135)  60  0.529412  473 

(255,133)  61  0.521569  481 

(255,131)  62  0.513725  489 

(255,129)  63  0.505882  497 

(255,127)  64  0.498039  505 

(255,125)  65  0.490196  513 

(255,123)  66  0.482353  521 

(255,121)  67  0.474510  529 


0.562500 

0.708333  -  -  T" 

0.781250  9  0.562500 

0.825000  9  0.450000 

0.854167  17  0.708333  9 

0.875000  17  0.607143  9™ 

0.890625  25  0.781250  9 

0.902778  25  0.694444  17 

0.912500  33  0.825000  17 

0.920455  33  "0J50000  TT 

0.927083  41  0.854167  25 

0.932692  41  0.788462  25 

0.937500  49  0.875000  25 

0.941667  49  0.816667  33 

0.945313  57  0.890625  33 

0.948529  57  0.838235  33 

0.951389  65  0.902778  41 

0.953947  65  0.855263  41 

0.956250  73  0.912500  41 

0.958333  73  0.869048  49 

0.960227  81  0.920455  ~ 

0.961957  81  0.880435  49 

0.963542  89  0.927083  57 

0.965000  89  0.890000  57 

0.966346  97  0.932692  57 

0.967593  97  0.898148  65“ 

0.968750  105  0.937500  65 

0.969828  105  0.905172  65* 

0.970833  113  0.941667  73 

0.971774  113  0.911290  73" 

0.972656  121  0.945313  73 

0.973485  121  0.916667  81 

0.974265  129  0.948529  81 

0.975000  129  0.921429  81 

0.975694  137  0.951389  89 

0.976351  137  0.925676  89 

0.976974  145  0.953947  89 

0.977564  145  0.929487  ~ 

0.978125  153  0.956250  97 

0.978659  153  0.932927  97 

0.979167  161  0.958333  105 

0.979651  161  0.936047  10? 

0.980114  169  0.960227  105 

0.980556  169  0.938889  113 

0.980978  177  0.961957  113 

0.981383  177  0.941489  113 

0.981771  185  0.963542  121 

0.982143  185  0.943878  121 

0.982500  193  0.965000  121 

0.982843  193  0.946078  129 

0.983173  201  0.966346  129 

0.983491  201  0.948113  12? 

0.983796  209  0.967593  137 

0.984091  209  0.950000  757 

0.984375  217  0.968750  137 

0.984649  217  0.951754  145 

0.984914  225  0.969828  145 

0.985169  225  0.953390  145 

0.985417  233  0.970833  153 

0.985656  233  0.954918  153 

0.985887  241  0.971774  153 

0.986111  241  0.956349  161 

0.986328  249  0.972656  161 

0.986538  249  0.957692  161 

0.986742  257  0.973485  169 

0.986940  257  0.958955  16? 


0.562500 

0.482143 

0.421875 

0.708333 

0.637500 

0.579545 

0.781250 

0.721154 

0.669643 

0.825000 

0.773438 

0.727941 

0.854167 

0.809211 

0.768750 

0.875000 

0.835227 

0.798913 

0.890625 

0.855000 

0.822115 

0.902778 

0.870536 

0.840517 

0.912500 

0.883065 

0.855469 

0.920455 

0.893382 

0.867857 

0.927083 

0.902027 

0.878289 

0.932692 

0.909375 

0.887195 

0.937500 

0.915698 

0.894886 

0.941667 

0.921196 

0.901596 

0.945313 

0.926020 

0.907500 

0.948529 

0.930288 

0.912736 

0.951389 

0.934091 

0.917411 

0.953947 

0.937500 

0.921610 

0.956250 

0.940574 

0.925403 

0.958333 

0.943359 

0.928846 

0.960227 

0.945896 


(255,117)1  69  1  0.458824|  545  |  0.987319|  265  0.960145  177  0.961957 

(255,115)  70  0.450980  553  0.987500  273  0  975000  177  0.948214 

(255,113)  71  0.443137  561  0,987676  273  0.961268  177  0.934859 

(255,111)  72  0.435294  569  0.987847  281  0.975694  185  0.963542 

(255,109)  73  *0.427451  ~  577  0.988014  281  0.962329  185  0.950342 

(255,107)  74  0.419608  585  0.988176  289  0  976351  185  0,937500 

(255,105)"  ~  75  ‘0411765  593  ~  0.988333  289  0.963333  193  0.965000 

(255,103)  76  0.403922  601  0.988487  297  0.976974  193  0.952303 

(255,101)' ~~  77  ~a396078  609  ~  0.988636  297  0.964286  193  0.939935 

(255,99)  ~  78 '  ~  0.388235  617  0.988782  305  0.977564  201  0.966346 

(255,97)  ~  79 .  0.380392  625  0.988924  305  0.965190  201  0.954114 

(255,95)  80  0.372549  633  0.989063  313  0.978125  201  0.942188 

(255,93)  81  0.364706  641  0.989198  313  0.966049  209  0.967593 

(255,91)  82  0.356863  649  0.989329  321  0.978659  209  0.955793 

(255,89)  83  0.349020  ~  657  0.989458  321  0.966867  209  0.944277 

(255,87)  84  0.341176  665  0.989583  329  0.979167  217  0,968750 

(255,85)  85  0.333333  673  0,989706  329  0.967647  217  0.957353 

(255,83)  86  0.325490  681  0.989826  337  0.979651  217  0.946221 

(255,81)  87  0.317647  689  0.989943  337  0.968391  225  0.969828 

(255,79)  88  0.309804  697  0,990057  345  0.980114  225  0.958807 

(255,77)  89  0.301961  705  0.990169  345  0,969101  225  0.948034 

(255,75)  90  0.294118  713  0.990278  353  0.980556  233  0.970833 

(255,73)  91  0.286275  "  721  0.990385  353  0.969780  233  0.960165 

(255,71)  92  0.278431  729  0.990489  361  0.980978  233  0.949728 

(255,69)  93  0.270588  737  0.990591  361  0.970430  241  0.971774 

(255,67)  ~  94  *0.262745  745  0.990691  369  0.981383  241  0.961436 

(255,65)  95  0.254902  753  0.990789  369  0.971053  241  0.951316 

(255,63)  96  0.247059  761  0.990885  377  0.981771  249  0.972656 

(255^61)  ~  97*  ~o!2392TF  769  0.990979  377  0.971649  249  0.962629 

(255,59)  98  '0231373  777  0.991071  385  0.982143  249  0.952806 

(255,57)"  ~  99~  0.223529  785  0.991162  385  0.972222  257  0.973485 

(255,55)  100  0,215686  793  0.991250  393  0.982500  257  0.963750 

(255,53)  Tol  0.207843  801  0.991337  393  0.972772  257  0.954208 

(255,51)  102  0.200000  809  0.991422  401  0,982843  265  0,974265 

(255,49)  103  *0*792157  817  0.991505  401  0.973301  265  0.964806 

(255,47)  104  0.184314  825  0.991587  409  0.983173  265  0,955529 

(255,45)  105  0.176471  833  0.991667  409  0,973810  273  0.975000 

(255,43)"' ~  106  7)168627  841  0.991745  417  0.983491  273  0.965802 

(255,41)  107  0.160784  849  0.991822  417  0,974299  273  0.956776 

(255|39)"  ""  108™ '"0752941"  857  0.991898  425  0.983796  281  0.975694 

(255,37)  '  109  0.145098  865  0.991972  425  0.974771  281  0.966743 

(255,35)  TTo  0.137255  873  0.992045  433  0.984091  281  0.957955 

(255,33)  111  *0.129412  881  0.992117  433  0.975225  289  0.976351 

(255,31)  U2  0.121569  889  0.992188  441  0.984375  289  0.967634 

(255^29)  113  “0113725  897  0.992257  441  0.975664  289  0.959071 

(255,27)  114  0.105882  905  0.992325  449  0,984649  297  0.976974 

(255,25)  115  0.098039  913  0,992391  449  0,976087  297  0.968478 

(255,23)  116  0.090196  '  921  0.992457  457  0.984914  297  0.960129 

(255,21)  117  0.082353  929  0.992521  457  0,976496  305  0.977564 

(255,19)  118  0.074510  937  0.992585  465  0.985169  305  0.969280 

(255,17)  119  0.066667  945  0.992647  465  0.976891  305  0.961134 

(255,15)  120  0.058824  ~  953  0.992708  473  0.985417  313  0.978125 

(255^13)  ~  121  0050980  961  0.992769  473  0.977273  313  0,970041 

(255,11)  122  0.043137  969  0.992828  481  0.985656  313  0.962090 

(255,9)  123  0.035294  977  |  0.992886 1  481  0.977642  321  0.978659 

(255,7)  124  0.027451  985  0.992944  489  0.985887  321  0,970766 

(255,5)  125  0.019608  993  0.993000  489  0.978000  321  0,963000 

(255,3)  126  0.011765  1001  0.993056  497  0.986111  329  0.979167 

(255,1)  127  0.003922  1009  0.993110  497  0.978346  329  0.971457 
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6.  Conclusions  And  Plans  For  Future  Activity 


Program  participants  have  remarked  frequently  on  the  success  of  the  RAPCEval  contract  and 
on  the  quality  of  the  results  that  have  been  produced  to  date.  The  cooperation  among  various 
engineering  communities  and  their  fruitful  interaction  has  been  a  very  satisfying  byproduct  of  the  entire 
project.  Students,  university  faculty,  government  and  private  sector  engineers  have  united  in  the 
common  attack  on  a  number  of  priority  Air  Force  electronic  warfare  concerns.  The  student  master's 
degree  research  has  focused  on  topics  that  have  immediate  impact  on  the  transitioning  of  new  and 
improved  software  and  hardware  technologies  into  fielded  systems. 

Recently,  several  students  have  had  the  privilege  of  presenting  their  work  at  national  and 
international  conferences,  where  comments  have  been  quite  favorable.  Ron  Brinkley  presented  a  paper 
entitled  ‘Reed-Solo man  Coding  and  GMSK  Modulation  for  Cellular  Digital  Packet  Systems”  at  an 
IEEE  Pacific  Rim  Conference  on  Communications,  Computers,  and  Signal  Processing  (PACRIM’  97). 
Henderson  Benjamin  gave  a  talk  entitled  “Selection  of  the  Most  ‘Efficient’  Reed-Solomon  Code  for  a 
Specific  Application  Using  Neural  Networks”  at  the  International  Conference  on  Communications 
(ICC  ’99),  in  Vancouver,  British  Columbia.  These  accomplishments  are  appreciated  all  the  more  when 
it  has  been  noted  that  it  is  very  unusual  for  a  master’s  level  graduate  student  to  present  a  paper  at 
national  and  international  conferences. 

Twelve  master's  degrees  have  been  awarded  at  the  time  of  this  report.  The  title  of  each  report 
is  listed  in  the  table  together  with  availability  information  for  each  document. 


Author 

Report  Title 

Availability 

Mark  Astin 

“Application  of  Parallel 
Computing  Techniques  to 
the  RAD  Algorithm” 

Classified 

available  from  MERC  or 

AFRL/SNRP: 

#AFRL-SN-WP-TR- 1 998- 1 088 

Henderson 

Benjamin 

“Neural  Network  System 
the  Selects  Reed-Solomon 
Codes  for  a  Specific 
Application” 

Unclassified 

available  from  MERC  or 

AFRL/SNRP: 

#AFRL-SN-WP-TR- 1 998- 1056, 

Section  5.2.9 

Ron  Brinkley 

“Burst  Error  Correction 
with  Reed-Solomon  Codes” 

Unclassified 

Submitted  to  AFRL/SNRP  for 
publication  -  anticipate  availability 

12/99 

Mark 

Campbell 

“Auto-Regressive  Spectral 
Analysis  -  EW 

Applications” 

Unclassified 

available  from  MERC  or 

AFRL/SNRP: 

#WL-TR-94- 1 057,  Apx.  E 
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Randy  Ford 

“Comparison  of 

Differential  Evolution  to 
the  Simplex  Method  in 
Optimization  During 

Passive  Emitter  Location” 

Unclassified 

Available  from  MERC 

In  process  of  submission  for 
publication  to  AFRL/SNRP 

Claus 

Franzkowiak 

‘Tour-Pulse  Primary  RAD 
Filter  Development” 

Classified 

available  from  MERC  or 

AFRL/SNRP: 

#AFRL-SN-WP-TR- 1 998- 1 087 

Neal  Gamer 

“Error  Correction  and 
Prediction  for  Improved 
Communication  of  Time 
and  Time  Measurements” 

Unclassified 

available  from  MERC  or 

AFRL/SNRP: 

#WL-TR-96- 1161,  Apx.  D 

Joseph 

Kelley 

“A  Parameter 

Determination  Alternative 
for  RAD  Analysis” 

Classified 

Available  from  AFRL/SNRP 
#WL-TR-95-1005 

Joseph 

Kelley 

“MultiGroup  Simultaneous 
RAD  Parameter  Selection” 

Classified 

Available  from  AFRL/SNRP 
#WL-TR-97- 1 094 

Max  Roesel 

“Agile  RF/PRI  Radar 
Analysis  via  RAD” 

Classified 

Available  from  AFRL/SNRP 
#WL-TR-95-l  020 

Dave  Schuler 

“Comparison  of  Algorithms 
for  Geolocation  of  Radar 
Signals” 

Unclassified 

Available  from  MERC- 
requires  need-to-know 

Kirk  Wright 

“Object  Oriented  Modeling 
of  the  AN/ALQ-172” 

Classified 

available  from  MERC  or 

AFRL/SNRP: 

#AFRL-SN-WP-TR- 1 998- 1086 

At  the  present  time,  several  students  are  planning  to  start  new  research-  Reports  from  various 
colleagues  in  the  military  and  industry  indicate  that  additional  students  will  be  interested  in  participating 
in  projects  of  the  kind  sponsored  by  the  RAPCEval  program. 
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